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ABSTRACT 

We present a new maximum likelihood method for the calculation of galaxy luminos- 
ity functions from multi-band photometric surveys without spectroscopic data. The 
method evaluates the likelihood of a trial luminosity function by directly comparing 
the predicted distribution of fluxes in a multi-dimensional photometric space to the ob- 
servations, and thus does not require the intermediate step of calculating photometric 
redshifts. We apply this algorithm to ~ 27, 000 galaxies with ^ 25 in the MUSYC- 
ECDFS field, with a focus on recovering the luminosity function of field galaxies at 
z < 1.2. Our deepest LFs reach M r —14 and show that the field galaxy LF deviates 
from a Schechter function, exhibiting a steep upturn at intermediate magnitudes that 
is due to galaxies of late spectral types. 

Key words: galaxies: luminosity function, galaxies: evolution, methods: statistical 



1 INTRODUCTION 

The galaxy luminosity function (LF) has long been one of 
the cornerstones of statistical observational cosmology. As 
the distribution of galaxies as a function of luminosity is 
closely related both to the halo mass function predicted by 
structure formation models and to the astrophysical mech- 
anisms governing star formation, the LF provides a funda- 
mental test for models of galaxy formation and evolution. 
The low-redshift LF has also b e en ext ensively studied (see 
BingR'cli, Sand age fc TammanrJ |l988) for a review), and 
some convergence has been reached on its genera l shape at 
the br ight end. In the famous parametrization bv lSchechterl 
l|l976t ). a faint-end power-law slope of a ~ —1.2 is the con- 
sensus value. The characteristic bright-end absolute magni- 
tude M* is dependent on the magnitude system an d filter 
used, but the determination bv lBlanton et all ()2003h sets it 
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at M* = -21.18 ±0.04 in the r-band at z=0.1 (adjusted for 
#o = 71). 

Nevertheless, the LF remains a field of ongoing study 
as modern instrumentation allows us to survey increasingly 
fainter galaxies. One aspect that has received great atten- 
tion is the evolution of the LF with redshift, as this provides 
a much more stringent test of galaxy evolution models than 
the single data point that the local LF offers us. Several 

1 (|Wolf et al 
3 (iPoli et al 



group s have studied the LF to limits of z f 
2003) or even higher, up to and beyond z 
20031 : iGiallongo et alj 120051 : iMarchesini et al. I 120071 ). and 



have found evidence for evolution. Broadly speaking, early- 
type galaxies show a marked number density evolution with 
redshift, while late- type galaxies exhibit a luminosity evolu- 
tion. 

But another frontier in our exploration of the LF 
lies at low redshifts and very faint absolute magnitudes. 
Traditionally, the deepest probes of the LF come from 
studies of galaxy clusters. Providing a sample of hundreds 
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of galaxies at the same distances, clusters offer an op- 
portunity to determine LFs simply by counting galaxies 
and then applying corrections for sample contamination 
by unrelated background galaxies. This meth od has been 
criticized by IValotto. Moore fc Lambasl l|200l[ ) and shown 
to potentially produce artificial ly steepened faint-en d slopes 
for optically-selected clusters (|Valotto et alJ l2004h due to 
the fact that optical selection favors clusters with a higher- 
than-average density of background galaxies. However, 
since it eliminates the need for spectroscopic observations 
to determine the exact distance of each galaxy, this method 
has nonetheless been widely applied and has probed the 
LF of nearby clusters to spectacular depths ((Trenthaml 
1 19981 ; iTrentham fc HodgkinlEffloj iTrentham fc Tulry||2002l )T 
Frequently, LFs calculated in this way show an upturn in 
the LF, setting in several magnitudes below the character- 
istic bright-end magnitude M*, from a slope of a w —1.2 
to a much steeper slope of a ~ —1.5. M ore recently, a 
simila r upturn has also been claimed by IPopesso et al.l 
(2005) based on applying the background subtraction 
technique to clusters in the SDSS. However, for a sample of 
clusters at z ~ 0.3, th i s has also recently been debated by 
lHarsono fc de Proprisl ((2008), who are also working with 
statistical background subtraction, but find no evidence for 
an up turn in the cluster LF. Furthermore, iRines fc Gelleri 
(2008), working with a spectroscopic sample, find that the 
LFs of the Virgo Cluster and Abell 2199 can be represented 
by single Schechter functions consistent with the canonical 
values of a in the range of —1.1 to —1.3, and attribute 
the discrepancy from other studies to the use of statistical 
background subtraction without spectroscopic membership 
confirmation in the earlier work. In principle, an upturn 
is not implausible. It is known that galaxies with different 
morphologies or spectral types have very different Schechter 
functions; for example, the LF of late-type galaxies is 
generally very steep, while that of early-type galaxies is 
shallower. If both are extrapolated to faint magnitudes 
and then co- added, a faint-en d upturn can be reproduced 
qualitatively ijWolf et al.ll2003h . provided that the normal- 
ization of the late-type LF is still high enough in these 
dense environments in order for it to dominate the overall 
LF at some point. Some surveys of nearby groups and 
clusters have also repor ted a "dip" in the LF at interme- 
diate mag nitudes (|Flint . Bolte fc Mendes de Qliveiral 

20031: iMiles. Ravchaudhurv fc Russell l200rj; 

Mendes de Oliveira. Cypriano fc Sodrel 1200a : iBiviano et al.l 



1996), which could also be an indication of a bimodality 



that would be revealed as the superposition of a shallow 
and a very steep Schechter function in a sufficiently deep 
and large survey. However, the fact that spectroscopic 
surveys of the field or clusters have not yet unambiguously 
verified this result (because of their generally shallower 
depth) makes this observation controversial. 

It would be extremely important to know if this up- 
turn is real and if it is specific to cluster environments (and 
thus an indicator of environment-specific processes generat- 
ing large numbers of dwarf galaxies) or rather a universal 
feature of the LF. Current Cold-Dark-Matter cosmologies 
surmise that galaxies are embedded in dark matter halos, 
and the low-mass end of the halo mass function is thought 
to be quite steep, with a power-law slope of a m —1.8. 
This is in marked contrast with the LF of field galaxies 



in the local unive rse, which exhibits a fairly flat faint-en d 
slope of a « -1.2 jBlanton et al.ll200ll ; lNorberg et al]|2002l ). 
However, the field LF is well-established from spectroscopic 
surveys only to intermediate magnitudes of M T ~ —17. 
Also, studies of the local LF in the only environment where 
we can probe its extreme faint end directly, the Local 
Group and the nearest of galaxy groups beyond it, pre- 
fer a flat faint end slope a ll the way down to the extreme 
dwarf regime of M r m — 10 dPritchet fc van den Berghlll999l : 



iFlint. Bolte fc Mendes de Oliveirall2003h or suggest at most 



a mo derately steep slope I Trentham. Sampson fc Baneriil 
120051 ): these studies, however, suffer from being able to ac- 
cess only a very small cosmic volume. Although new Local 
Group members conti nue to be found up to the present day 
l|Simon fc G cha 2007), it is unlikely that they will change 
the census substantially at any but the very faintest lumi- 
nosities, substantially fainter than the limit at which the 
steep faint-end upturn in clusters has been claimed to ap- 
pear. Mechanisms su ppressing the format ion of galaxies in 
such low-mass halos l|Koposov et al.l (20091 ) may explain the 
small number of observable dwarf galaxies at low luminosi- 
ties, but it is unclear how such mechanisms would account 
for the presence of a large dwarf galaxy population in clus- 
ters. 

A tentative suggestion of an upturn in the fa int end of 
the field LF was found by iBlanton et al.l l|2005l ). who de- 
tected a steepening of the slope to a = —1.3 at intermedi- 
ate luminosities, but speculated that inclusion of low-surface 
brightness galaxies undetectable by the SDSS might modify 
this value toa~ —1.5 or even steeper. Apart from this work, 
the aforementioned studies of the LF in clusters provide the 
only suggestion of a steeper faint-end slope so far, but these 
are questionable due to their samples lacking spectroscopic 
confirmation and therefore being susceptible to background 
contamination. 

The frontier in exploring the LF at higher redshift 
and to fainter absolute magnitudes therefore lies in push- 
ing the flux sensitivity of samples towards fainter and 
more distant objects. However, at such low fluxes, spec- 
troscopic redshift determination remains extremely time- 
consuming and has a high failure rate, precluding its appli- 
cation to large samples. As an alternative approach, pho- 
tometric redshifts are increasingly being substit uted for 
spect r oscopic redshifts even in the study of LFs l|Benitezl 
l200d : iBolzonella. Pello fc Maccagnil |2002j ). However, it is 
well-known that photometric redshifts are not an equiva- 
lent substitute for spectroscopic redshifts; their accuracy is 
typically of the order of Az = 0.05(1 + z) even in the best 
cases, and frequently worse than that, leading to the galaxy 
distribution being smeared out in luminosity and distance. 
In addition, they suffer from catastrophic failures; galaxies 
of very different spectral types and at different redshifts can 
have similar colors, leading to large errors in their recov- 
ered photometric redshift that do not obey Gaussian errors. 
Even purely Gaussian errors in the photometry therefore 
inevitably propagate into non-Gaussian errors in the pho- 
tometric redshift, and samples selected by photometric red- 
shifts are inevitably contaminated by galaxies from galaxy 
populations outside these selection limits. In addition, the 
distribution of these errors varies from galaxy to galaxy, de- 
pending on its spectral properties, true redshift, and magni- 
tude. For all these reasons, the recovered photometric red- 
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shifts and, by extension, LFs, are always dependent on (ex- 
plicit or implicit) prior assumptions about the LFs of all 
galaxy populations at all redshifts — assumptions which, 
except in the case of very careful analyses, are often ne- 
glected, or in conflict with the recovered results. 

Unfortunately, luminosity function algorithms have 
adapted only poorly to this situation. The algorithms ap- 
plied to the problem of recovering LFs from photometric- 
only samples are based on the same methods used twenty 
years ago for the recovery of LFs from samples with known 
spectroscopic redshifts. In the simplest procedure, photo- 
metric redshifts are simply substituted for spectroscopic 
redshifts, and, assuming that the errors in the photo- 
metric redshifts are sufficiently understood, Monte Carlo 
simulations can be used to estimate whether a serious 
bias has been introduced into the recovered LF parame- 
ters. This approach can be successful, especially at high 
redshift, where a given uncertainty Az in the photomet- 
ric redshift translates into only a relatively minor uncer- 
tainty in the absolute magnitude. An example of a suc- 
cessful application of this procedure to t he Multiwave- 
length Survey by Yale-Chile ("MUSYC") (iGawiser et all 
20061) l|http://www.astro.yale.edu/MUSYC[ ) is found in 
Marchesini et all l|2007l ). Howe ver, at low redsh ift, as has 
been demonstrated in the past (|Chen et al.ll2003h . the prob- 
lems associated with the redshift uncertainties are generally 
larger: a redshift uncertainty Az ~ 0.1, typical for photo- 
metric redshift determinations, translates into a significant 
uncertainty in absolute magnitude. Correspondingly high- 
quality multi-ban d photometric da ta is needed to constrain 
the redshifts well (|Wolf et al.ll2003l ), but obtaining such ex- 
tensive datasets is not always practical. 

There have been suggestions in the literature for ad- 
justments to the LF calculation that take into account the 
photometric uncertainties and the resulting uncertainties 
in the re dshift. These have either emp l oyed Monte Carlo 
methods (|Bolzonella. Pello fc Maccagnil 120021 ) or incorpo- 
rated the photo-z uncertainties analytically into the proce- 
dure used for recovering the LF (|Chen et alj|2003h . However, 
both methods are based on the idea of evaluating a redshift 
probability distribution that is centered on the best-fit re- 
covered photometric redshift, and not centered on the true 
redshift. The net effect is that the galaxy distribution is sub- 
jected to a further convolution with the redshift error func- 
tion, instead of a deconvolution. 

The idea behind the Photometric Maximum Likelihood 
method (hereinafter: PML), which we introduce in this pa- 
per, is therefore to abandon the concept of assigning indi- 
vidual luminosities and redshifts to individual galaxies, since 
neither are known in non-spectroscopic surveys. Instead, we 
constrain the LF by requiring it to reproduce the distribu- 
tion of galaxies in a parameter space consisting of quantities 
that are known exactly: the observed fluxes in multiple filter 
bands. We refer to this parameter space, wherein observer- 
frame fluxes in multiple filter bands are used as the coordi- 
nates, as photometric space. This algorithm achieves a true 
deconvolution of the observed galaxy distribution into its 
constituent galaxy luminosity functions, and it allows for 
a self-consistent solution by avoiding priors that the recov- 
ered photometric redshifts and luminosity distributions in 
photometric surveys normally depend on. 

In this paper, we present the first application of the 



PML to a data set that is an ideal test case for this algo- 
rithm, the MUSYC survey (jGawiser et al.ll2006l ). MUSYC 
is a deep square degree survey (divided into four separate 
fields), comprising some ~ 200,000 galaxies. Photometric 
data in a multitude of optical and NIR filter bands exist. For 
the present study, we use one of the four MUSYC fields, the 
extended Chandra Deep Field South (E-CDFS), observed 
in eight filter bands: U, B, V, R, I, z, J, and K. The large 
number of filter bands make this survey very suitable to 
photometric redshift techniques as well as to the PML. The 
limiting depth of R > 26 mag means that we can probe the 
LF to extremely faint limits; the relatively wide area, com- 
pared to other pencil-beam surveys, guarantees that even 
the bright end will still be fairly well-sampled at all but the 
smallest redshifts. With a sample cut at R = 25, the PML 
permits us to constrain the field galaxy LF to a limit as faint 
as M r = -14. 

This paper is organized as follows: In jj^l we briefly dis- 
cuss the data set, the MUSYC survey. In fj3] we present our 
LF algorithm the Photometric Maximum Likelihood method 
and test it using a mock catalogue. In we present an 
application of our method to MUSYC by calculating lumi- 
nosity function parameters, split up by spectral type, from 
z « 0.05 to z ^ 1.2. 



2 DATA 

MUSYC is designed to provide a fair sample of the universe 
for the study of the formation and evolution of galaxies and 
their central black holes. The core of the survey is a deep 
imaging campaign in optical and near-infrared passbands of 
four carefully selected 30'x30' fields. MUSYC is one of few 
surveys to offer its combination of depth and total area, for 
additional coverage at X-ray, UV-, mid-infrared, and far- 
infrared wavelengths, and for providing the UBVRIz 1 JHK 
photometry needed for high-quality photometric redshifts 
over a square degree of sky (availability of NIR data varies 
from field to field). The primary goal is to study the prop- 
erties and interrelations of galaxies at a single epoch cor- 
responding to redshift ~ 3, but the favorable combination 
of depth, area, and passband coverage makes it suitable for 
studies of the general galaxy population over a wide range of 
redsh ifts as well as for Galactic astronomy l|Altmann et al.l 
120051 ). 

Imaging data for three of the four MUSYC fields was 
obtained with MOSAIC II on the 4 m Blanco Telescope 
at the Cerro Tololo Inter-American Observatory (CTIO) 
in Chile. For the fourth field, the Extended Chandra 
Deep Field South (E-CDFS), which we analyze in the 
present paper, our UBVRL imaging results from combin- 
ing public images taken with the ESO 2.2 m and WFI 
by the ESQ Deep Public Survey and COMBO-17 teams 
llErben et al. 1 120051: iHildebrandt et al. I l2005l : lArnouts et al.l 
l200ll ; IWolf et al.ll2004l ). Our z' imaging was taken with the 
CTIO 4 m and MOSAIC II. Our JK images of the E-CDFS 
were obtained with the CTIO 4 m and ISPI. 

Photometric coverage of the four MUSYC fields is not 
perfectly homogeneous. For the present study, we have thus 
chosen to focus on one of the four MUSYC fields, namely 
the E-CDFS, which offers the deepest optical imaging; fur- 
thermore, in contrast to the other MUSYC fields at the 
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time of this writing, the optically selected photometric cata- 
logue with UBVRIz photometry has been matched to NIR 
catalogues in J and K. The large number of spectroscopic 
redshifts available in the literature also enable us to check 
the quality of our models by estimating photometric red- 
shifts. This field is centered on RA = 03 : 32 : 29.0, 
Dec = -27 : 48 : 47 (J2000). 

The data that we use for this analysis are drawn from 
a catalogue selected by the combined B-, V-, and _R-band 
flux, and complemented by photometry in the J- and K- 
bands. The catalogues of J- and if-band photometric data 
for the BUi?-selected detections have not been published, 
but separate, isT-band-selected catalogues of this field based 
on the same imaging material c omprising the NIR b ands are 
already public and described in lTavlor et al.l (|2008h . 

The full catalog of the MUSYC E-CDFS contains 
84410 objects. At the depth of MUSYC, the detections 
consist primarily of galaxies, but the catalog also con- 
tains point sources (stars and QSOs), which we remove 
based on morphological criteria, using the Star/Galaxy clas- 
sifier (S/G) provided b y the Source Extractor software 
(Bcrt in fc Arnouti[l996), which was used in the construc- 
tion of the detection catalog. To calibrate the threshold 
value of S/G, we compare its distribution to the position 
of a detection in the thr-^ir — m,R parameter space. In 
such a parameter space, which is shown in Fig. [1] point 
sources populate an approximately horizontal band. Exper- 
imentally, we find that S/G = 0.97 separates the two pop- 
ulations with minimal overlap. This is demonstrated in Fig. 
Q]by marking all objects with S/G > 0.97 with asterisks in 
the left panel, and all other objects with dots in the right 
panel. This is a rather high threshold, compared to values 
adopted in the liter ature in simila r efforts (i.e., S/G = 0.9 , 
llCapak et al.ll2003h; S/G = .7, l|Rowan-Robinsonl |2004 ) ; 
S/G = 0.5. |Grutzbauch1l2005l )). and it is possible that our 
S/G threshold fails to exclude scattered light from saturated 
stars or close double star systems. To estimate the impact 
of the star/galaxy threshold, we have therefore carried out 
a separate analysis adopting a more conservative threshold 
of S/G — 0.8 and will compare its results with the default 
sample where appropriate in i|4] 

For the calculation of the LF from a flux-limited sample, 
a thorough understanding of the selection function is imper- 
ative. The optical MUSYC source catalogs are detected by 
SExtractor i n the combined B-, V-, and J?-band image, as 
described in lGawiser et a l. (2006). The limiting magnitude 
for the E-CDFS field is BVR = 27.1 mag (in the AB magni- 
tude system, which is used throughout this paper). For the 
present paper, we apply an additional cut in just a single 
filter band, the _R-band, to create our subsample. At R ^ 25 
mag, corresponding to a flux limit of R — 0.363 piJy, which is 
the limit adopted for our analysis, the sample is complete in 
BVR. Objects with SExtractor warning flags 1, 2 (presence 
of a close neighbour object possibly affecting the photome- 
try; necessity of deblending) are left in the catalog; objects 
with higher warning flags are removed without statistical 
compensation. 

For our analysis, we u se "corrected aperture magni- 
tudes" l|Gawiser et al.ll2006l ). which yield an estimate of the 
total flux, but offer good S/N for compact sources. An ex- 
ception are very extended objects (with a half-light radius 
in the BV R- image of 1 arcsec or more), for which we sub- 
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Figure 1. Difference between average surface brightness (in mag 
arcsec -2 ) and magnitude vtir versus rriR. Point sources tend 
to populate a horizontal branch in this diagram, and extended 
sources an extended cloud on the right side, but the two over- 
lap for faint magnitudes. In the left panel, objects with the 
star/galaxy class S/G > 0.97 are marked, in the right panel ob- 
jects with S/G < = 0.97. This limit decomposes the two popu- 
lations with minimal overlap; we therefore accept objects with 
S/G <= 0.97 as galaxies. 



stitute the "AUTO" photometric measurement provided by 
SExtractor. In addition, we have carried out modifications 
to the published MUSYC photometry, as described in Ap- 
pendix [X] to adjust the photometric zero points and ensure 
that the photometric errors conform to a Gaussian distri- 
bution. These corrections have significantly improved the 
quality of the recovered photometric redshifts in compari- 
son to a subsample of galaxies with redshift estimates from 
the literature, spanning a broad range of redshifts and SED 
types, and have brought the uncertainty distribution in line 
with the assumptions made in our algorithm. 



3 METHOD 

3.1 Principles of the Photometric Maximum 
Likelihood Approach 

The most common type of LF algorithms for magnitude- 
limited samples of field galaxies are maximum likelihood 
algorithms. These algorithms are based on the principle of 
calculating a predicted distribution of galaxies from a trial 
LF, comparing it to the observed distribution, determining 
the likelihood of the trial LF, and then modifying the latter 
until the likelihood has been maximized. Traditionally, the 
distributions that are being compared are the distribu- 
tions over luminosity llEfstathiou. Ellis fc Peterson! 1 19871 : 
ISandage. Tamman fc Yahillll979r ). sometimes over luminos- 
ity plus additional para meters that the LF could depend 
on, such as spectral type (jBromlev et al . 1998; Cr oton et al.l 
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20051 ') , morpholog y ( Binggeli. Sandage fc Tammannl 
198a i, or redshift l|Wolf et all 120031 ; iPoli et all 120031 ; 



Giallongo et~ai1l2005h itself. 



In photometric surveys, neither spectral type, nor red- 
shift, nor luminosity are known exactly. Nonetheless, past 
attempts to calculate the galaxy LF from photometric sur- 
veys have sought to reconcile the data with the algorithm by 
simply substituting photometric redshifts for spectroscopic 
redshifts. This can potentially cause several severe problems: 

• Random errors in the photometric redshifts "blur" the 
recovered distances and absolute magnitudes. This intro- 
duces systematic biases in the recovered LF, as can be il- 
lustrated by a gedankenexpenment: A typical LF is well- 
approximated by a power law at the faint end and has a 
well-defined exponential downturn at the bright end. If ran- 
dom errors exist in the galaxy redshift, the absolute magni- 
tude estimates of all sample galaxies will be convolved with 
an error function, which "blurs" the shape of the LF and 
moves the bright-end "knee" to brighter luminosities. This 
can be minimized by operating at high redshift (where the 
impact of redshift uncertainties on the absolute magnitude 
is smaller) and, of course, with high-quality multi-band pho- 
tometric data. 

• Systematic errors and catastrophic failures in the photo- 
metric redshifts may severely impact the recovered LF. For 
example, systematic offsets between the photometric red- 
shifts and the true redshift (a common occurrence resulting 
from errors in the photometric calibration) will directly in- 
troduce shifts in the recovered absolute magnitudes. Catas- 
trophic failures (i.e., photometric redshifts significantly and 
systematically different from the true redshifts) can have a 
devastating impact; for example, a small number of nearby, 
bright galaxies mistakenly fitted as high-redshift galaxies 
may easily dominate the bright end of the LF recovered at 
high redshift. It is important to note that such catastrophic 
failures can result simply from the propagation of Gaus- 
sian errors of the photometry into the photo-z determina- 
tion. This means that modeling photometric redshift uncer- 
tainties as Gaussians is not legitimate even for perfectly- 
calibrated photometry. 

• Sample contamination is common. LF determinations 
are usually carried out for subsamples of galaxies in a con- 
fined redshift range and for selected SED types only. How- 
ever, any subsample selected from photometric redshift data 
only will almost certainly be contaminated by galaxies at 
other redshifts and with a wide range of SED types, and 
incomplete with respect to galaxies that meet the selec- 
tion criteria (see for example iTavlor et al.l (|20Qg| ). which ar- 

ues that photom etric redshifts from the COMBO-17 survey 
Wolf et aH 120041 ) suffer from contamination at z > 1 due 
to the lack of NIR photometry). Aggravating this problem, 
the amount of contamination is, of course, dependent on the 
LF of the contaminant population. Therefore, even if some 
attempt is made to minimize the contamination (e.g., by 
introducing an LF prior into the photo-z calculation), the 
resulting LF is dependent on the LFs of all other galaxy 
populations at other redshifts. 

Various procedures have been suggested for propa- 
gating the inherent uncertainties from photometric red- 
shifts into the LF. A simple st rategy, proposed by 
iBolzonella. Pello fc Maccagnil i|2002l ), is to carry out a 



Monte Carlo simulation, drawing redshifts for each galaxy 
at random from the probability distribution indicated by the 
photo-z algorithm. However, the fallacy in this approach is 
that the Monte Carlo redshifts are drawn from a probability 
distribution that is centered on the best-fit photo-z, and not 
on the true redshift. While this approach does indeed propa- 
gate the photo-z uncertainties into the final solution, it does 
not correct the systematic biases in the photo-z recovery; 
on the contrary, drawing galaxies from a photo-z probabil- 
ity distributions around photo-zs that are already affected 
by these uncertainties will convolve the solution even more, 
instead of deconvolving it, and may even exacerbate the 
systematic e rrors. A slightly different approach is taken by 
IChen et al.l (|2003l ). who incorporate the photo-z probability 
distribution into the likelihood function itself. However, in- 
spection of their ansatz shows that the effect is, in fact, again 
a further convolution of the galaxy distribution, rather than 
a deconvolution, equivalent to the Monte-Carlo approach (it 
is equivalent to replacing each sample galaxy with a sum of 
"fractional" galaxies with a range of redshifts); it is thus 
a more elegant approach suffering from the same problem: 
the redshift probability distribution is evaluated around a 
redshift that is not the true redshift. 

A mathematically much more consistent method that 
addresses many of the criticism s above has been suggested 
more recently by IShethl (|2007f ) : The recovered photometric 
redshift is treated as a quasi-observable distinct from the 
true redshift, and reproducing its distribution is an addi- 
tional requirement imposed on the LF algorithm. This ap- 
proach appears theoretically sound; however, it stops short 
of replacing the photometric redshift with a different, and 
easier-to-handle, observable, as we propose here. This raises 
two practical problems for the application: 1) The algorithm 
is still based on the concept of a best-fit photo-z. The proce- 
dure of first calculating a best-fit photo-z and then feeding 
it as input into an LF algorithm always carries with it a 
loss of information that is contained in the full photomet- 
ric vector. 2) The error function of the photometric redshift 
must be modeled well, which is challenging, as it is a func- 
tion of many parameters, including the observed flux and 
true SED of a galaxy. The shape of the photometric red- 
shift probability distribution function varies with the SED, 
redshift, and apparent magnitude of the galaxy in question 
(e.g., late-type galaxies are much more difficult to constrain 
than early- type galaxies), and may not even generally be as- 
sumed to be Gaussian either, but may in some cases exhibit 
a complex shape with secondary maxima. 

An additional issue that has not been addressed by con- 
ventional LF al gorithms (but appears to be accounted for 
bv lShethl l|2007tl ) is that the probability distribution of pho- 
tometric redshift is itself dependent on implicit or explicit 
prior assumptions about the LF of galaxies at any redshift. 
Photometric redshifts are therefore dependent on prior as- 
sumptions about the redshift distribution of galaxies (i.e., 
the LF of galaxies of all possible spectral types throughout 
cosmic history). In photometric redshift determinations, this 
dependence is often neglected, and redshifts are calculated 
simply as best-fit estimates, regardless of the redshift or tem- 
plate type that yields the fit. Note, however, that even this 
approach to photo-z determination carries an implicit and 
unjustified prior, namely, that the prior probability distri- 
bution, which depends on the cosmic volume and the LF, is 
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flat with redshift. A m ore promising approach is Bayesian 
redshift determination (|Benitej|200C ). which applies priors 
(usually in the form of luminosity function priors) to select 
the most probable photometric redshift. However, while at 
least the shape of the bright end of the LF is reasonably 
well known, priors about the faint end and the high-redshift 
LF are highly uncertain. Furthermore, upon obtaining a new 
solution for the LF at any redshift, the priors should be mod- 
ified to accommodate this new information, and the process 
iterated, so as to avoid a contradiction between the LF result 
and the assumed prior. If this is not done, a logical inconsis- 
tency between the prior applied to the photo-z calculation 
and the resulting LF will generally result. 

In summary, there are four points of criticism that may 
be raised against the methods employed to recover LFs from 
photometric samples in the past: 

• Instead of a deconvolution, they have often carried out 
a further convolution of the data by applying an assumed 
error function to the observed data. A proper deconvolution 
requires that the error function be applied to a model, and 
the resulting distribution be compared to the observed data. 

• They have evaluated the predicted and observed galaxy 
distributions in a parameter space consisting of redshift and 
derived variables, which are not actual observables, but are 
affected by very complex error functions, which prohibit a 
direct deconvolution. In such cases, the comparison is more 
easily carried out in a parameter space consisting of direct 
observables, allowing us to convolve a predicted data set 
with the known error function instead. 

• They have used photometric redshifts and have conse- 
quently been forced to evaluate their error functions. Photo- 
metric redshifts have very complicated error functions that 
are usually neither Gaussian nor uniform, but instead fea- 
ture secondary maxima and are dependent on magnitude 
and SED type. But the calculation of LFs does not require 
knowledge of the redshift of individual galaxies. We there- 
fore suggest to eliminate the photometric redshift from the 
consideration and focus on direct observables, i.e., the ob- 
served fluxes, which have a much simpler error function. 

• They have often failed to take into account the de- 
pendence of the photometric redshift on prior assumptions 
about the LF at all redshifts. A proper determination of the 
LF must acknowledge this dependence and ideally provide 
a self-consistent solution that is not dependent on external 
priors, but constrained only by the data. 

The problems described above arise from the attempt 
to use established LF algorithms, which were originally de- 
veloped for samples for which the exact redshifts and, by 
extension, absolute magnitudes, are known, in a context 
that they were not constructed for, namely, photometric- 
only galaxy surveys. Our Photometric Maximum Likelihood 
method adopts the basic principle of maximum likelihood 
algorithms: From a trial LF, the galaxy distribution in a 
certain parameter space (e.g., absolute magnitude) is pre- 
dicted, and the probability for drawing the observed sample 
from this distribution (i.e., the likelihood of the trial LF) 
is then calculated. The trial LF is then adjusted to opti- 
mize this likelihood. However, the Photometric Maximum 
Likelihood algorithm differs from established algorithms by 
comparing the observed and predicted galaxy distributions 
in a parameter space consisting only of direct observables, 



namely, the observed photometric fluxes in all available fil- 
ter bands. We refer to this parameter space as photometric 
space. This addresses the problems described above. 

• The PML is able to carry out a proper deconvolution of 
the observed photometric data: It applies an error function 
to the expectation values for the photometric fluxes in each 
filter band and compares the resulting distribution to the 
observed data. 

• By evaluating the galaxy distribution as a function of 
observed fluxes, instead of photometric redshifts, the PML 
avoids the unnecessary loss of information that is usually in- 
curred when converting the full information of a multi-band 
photometric observation into a single photometric redshift. 
This also offers the decisive benefit that photometric errors 
in individual filter bands are much easier to model than 
those of photometric redshifts. 

• The algorithm considers the contributions of all LFs 
at all redshifts to the distribution in photometric space si- 
multaneously, and therefore finds a self-consistent solution 
for each LF that is (ideally) not dependent on implicit or 
explicit external luminosity priors. 

Our algorithm, by design, naturally accounts for the de- 
pendence of the LF on redshift and SED type. An important 
aspect is the fact that, when evaluating the distribution of 
galaxies in photometric space, in principle, any galaxy pop- 
ulation at any redshift can contribute to the galaxy counts 
at a given point in photometric space. It is not possible to 
clearly separate samples by SED type and redshift. This 
implies that the likelihood found by the PML algorithm is 
calculated as a likelihood for all LFs of all galaxy popula- 
tions throughout cosmic history. This appears at first as a 
drawback, greatly enlarging the number of free parameters 
involved in the problem. However, as has been pointed out 
above, in the conventional approach, it is not strictly possi- 
ble to clearly separate a sample by redshift and SED type 
either. Any association of an individual galaxy to a given 
redshift and SED type is only a statistical one, and under 
proper Bayesian consideration, the probability for any such 
association is dependent on the prior assumptions about the 
LFs of all galaxy populations that a galaxy could in princi- 
ple have been drawn from. Whereas this dependence usually 
goes unacknowledged, it is made explicit in the PML, and 
the problem of introducing external priors is (in an ideal 
case) eliminated by solving for all contributing LFs simul- 
taneously. 

3.2 Derivation 

The derivation of a maximum likelihood algorithm starts 
with the definition of the likelihood function, i.e., the prob- 
ability that the observed data of a galaxy survey would have 
been obtained from a given trial LF. The specific likeli- 
hood function t hat w e use is inspired by the approach of 
iMarshall etaD (Il983l). which has also been used more re- 
centlv by iPoli et all (|2003l ) and iGiallongo erall (|2005l ) for 
the determination of galaxy LFs, which uses a product of 
Poisson probabilities to derive the survey probability. 

We imagine an n-dimensional flux space, the photomet- 
ric space. In this space, each dimension corresponds to one 
filter band in which a flux is measured, and each coordinate 
to a measured flux in the corresponding filter. Each galaxy i 
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with measured photometry thus occupies a unique position 
fi that represents the observed flux in all n filters. We now 
imagine this space to be subdivided into small cells [f;df] 
that are chosen small enough that each cell contains either 
exactly or 1 galaxies. We then define the likelihood func- 
tion of the survey as the product of the Poisson probabilities 
that each cell that contains a galaxy would contain exactly 
one object, and that each of the remaining cells would con- 
tain none, under the prior of a chosen luminosity function. 
In other words, 



-Hfi)df 



-Hfj)dJ 



(1) 



Here, A(/;) df is the expectation value of the number 
of galaxies in cell [/; ;df] of photometric space. The index i 
runs over all observed galaxies, while the index j runs over 
all cells of fl ux space that are not occupied by a g alaxy. 
In contrast tolMarshall et alJ (|l983h . IPoli et al.1 (120031 ). and 
iGiallongo et al.l l|2005l ). we apply this likelihood evaluation 
to the distribution of galaxies in photometric space, while 
the previous work makes use of the availability of full or par- 
tial spectroscopic coverage to evaluate the galaxy distribu- 
tion as a function of luminosity and redshift. This distinction 
is the defining feature of the PML. 

Recently, Kelly, Fan fc Vestergaardl (|2008h have sug- 
gested that the proper statistical distribution to use is a 
binomial one, rather than the Poisson distribution used in 
the ansatz of Marshall et al. However, they acknowledge 
that, for surveys that sample only a small fraction of the 
total sky, which is certainly the case for MUSYC, Poisson 
distributions are acceptable. 

In our sample, we do not have the luxury of spectro- 
scopic redshifts. It is therefore necessary to introduce a map- 
ping between the luminosity function parameters and photo- 
metric space. We assume that the photometric data points 
are distributed according to a Gaussian distribution cen- 
tered on the expectation value for a given absolute magni- 
tude Mo, SED, and redshift zo, and furthermore, that the 
photometric errors are independent in each filter. 

We do this by defining A in the following way: 



Hh) = 



SED 



x<E>(M ; SED; z | P). (2) 



This expression integrates over all possible galaxy pop- 
ulations — at all redshifts, all absolute magnitudes, and all 
spectral types SED — and evaluates their contributions to 
the galaxy density in photometric space at the flux position 
fi. The term dV c /dz is the differential comoving volume, 
and P are the assumed parameters of the LF whose like- 
lihood we wish to test. The fraction of galaxies with the 
physical parameters Mo, SEDq, and zo that contribute to 
the galaxy density at a point /; is determined via the prob- 
ability term p a — the \ 2 probability density resulting from 
a comparison between the fluxes predicted for these param- 
eters, and the observed fluxes; this distribution is used here 
under the assumption that the errors in each filter band are 
uncorrelated with each other and that the photometric un- 
certainties in each band are Gaussian. This probability term 
is our assumed photometric error function. The denomina- 



tor is the product of the observed flux uncertainties and 
serves to normalize the probability distribution. We use the 
following definition of the dimensionless probability density: 



P* = 



(2tt) D / 2 



-X 2 /2 



(3) 



Here, D is the number of degrees of freedom (in our case, 
the number of filters with valid photometry). 
The x 2 quantifier is defined as 



D 

£ 

d 
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Id 



Ml 



(4) 



We use the empirical measurement uncertainties Af d of 
galaxy i to define the \ 2 statistic. This is a compromise. A 
more consistent approach would be to select uncertainties 
appropriate for the predicted fluxes of the galaxy in each 
filter band. However, this would require an error model re- 
producing accurate uncertainties in every filter band as a 
function of flux. As the photometric uncertainties are largely 
dictated by conditions other than luminosity alone (such as 
the surface area of the galaxy), this would require additional 
assumptions that would weaken, rather than strengthen, our 
case. We believe that the substitution of the empirical un- 
certainties in this definition is a justifiable compromise. 

A notable point about this approach is that the absolute 
magnitude of an individual galaxy enters via the normaliza- 
tion of the SED in any desired rest-frame filter band (which 
need not be identical to the filters used in the survey). The 
fc-correction is therefore integrated into the algorithm im- 
plicitly, and the absolute magnitude for a given galaxy at a 
given redshift and SED is constrained by all available filter 
bands, rather than by a single band and a corresponding k- 
correction; it is therefore usually quite precisely constrained. 
A caveat is that the accuracy of absolute magnitudes de- 
pends on the accuracy of the templates; therefore, although 
the algorithm is in principle capable of calculating LFs in 
any desired rest-frame filter band, it is desirable that the lu- 
minosity in this filter band be well-constrained by the avail- 
able photometry (i.e., we would distrust rest-frame LFs in 
filter bands at the extreme red or blue end of the survey 
coverage) . 

As Eq. [2] shows, the PML accounts intrinsically for 
the dependence of the LF on redshift, spectral type, and, 
if desired, other parameters of the LF. It also automati- 
cally accounts for the effect of photometric uncertainties 
on the recovered LF. Most importantly, however, it com- 
bines the functionality of a photometric redshift algorithm 
with that of a LF algorithm via the probability term Pa(fi \ 
Mo; SEDq; zo), which relates a given galaxy population to 
the galaxy distribution in photometric space. 

Equation [T] can now be transformed by setting S = 
—2 In C (and writing Xi as shorthand for A(/;)): 



S= -2 V] ZnAief e~ Ai<ff -2^2lne~ 



We split the argument of the first logarithm: 



S=-2V ln\ z df - 2 ^ Zr 



-Xidf 



(•») 



(6) 



8 D. Christlein et al. 



N N 

S = -2^ln\idf + 2^\idf + 2^2\jdf (7) 

i i j 

The second and third terms can be combined into an 
integral over all of observable flux space: 

N 

S = -2^2 lnX ^f + 2 Xdf (8) 

i 

Here, the first term represents the likelihood of drawing 
the observed galaxies from the assumed input LFs, while the 
second term is simply the total predicted number of galaxies 
within the flux limits of the survey. The expression df, the 
volume of one of our cells in photometric space, can be omit- 
ted, as it will only add a constant offset of S. Increasing the 
assumed space density of any galaxy population (for exam- 
ple, by increasing the normalization of the LF) will decrease 
the first term, reflecting the greater probability of drawing 
the observed galaxies from such an input LF, but it will also 
increase the second term, the total number of galaxies ex- 
pected in the survey. The luminosity function that balances 
these terms by minimizing the function S is the maximum 
likelihood LF. 

Calculating A(/) in the second term requires an as- 
sumption regarding the width of the distribution, i.e., the 
photometric errors of a hypothetical galaxy at a point / in 
photometric space, in order to correctly represent the con- 
tributions of galaxy populations near the flux limit of the 
sample to the total number of observed galaxies. It should 
be set to a value typical of the survey under consideration; 
we have experimentally found that the impact of this value 
on the recovered LFs is very minor. We set it to 25% of the 
flux, which is characteristic of most objects near this flux 
limit. This value determines the fraction of galaxies whose 
observed fluxes scatter across the survey limit, and it will 
influence the solution only if the faint end of the LF is very 
steep, i.e., if a much larger number of normally unobservable 
galaxies scatter above the survey limits due to photometric 
errors than vice versa. 

Despite our magnitude threshold of R ^ 25, there are 
galaxies in the catalogue that are consistent with zero flux in 
at least one other filter band. As a result, the algorithm may 
associate them with extremely faint populations that would 
normally not be expected to contribute to galaxy counts at 
all, given our model SEDs and the nominal flux error of 
25%. Inclusion of such objects may prevent a full conver- 
gence of the LF solution, especially in parameters that are 
constrained by few other objects. To prevent this, we only 
evaluate contributions to Ai from populations that yield a 
non-zero contribution to the total expectation value for the 
number of galaxies in the survey. In other words, we ig- 
nore potential contributions to the observed galaxy sample 
from populations that should, assuming reasonable SEDs 
and photometric errors, be outside the survey flux limits 
anyway. 

The fact that no precise redshift information is available 
does imply that the PML shares a disadvantage with the 
traditional V/Vmax method: It requires certain assumptions 
about spatial homogeneity of the LF, violation of which may 
bias the results. In a non-spectroscopic sample, this is gener- 
ally unavoidable, because photometric redshift uncertainties 
wipe out detailed information about large scale structure. To 



the extent that the photometric information permits it, the 
PML may constrain the LF in several independent redshift 
bins and thus reproduce gross features of the cosmic galaxy 
density fluctuation. 

3.3 Parametrization of the LF 

A number of possible parameterizations have been proposed 
in the literature to represent the galaxy luminosity func- 
tion. The classical representation is the Schechter function 
ISchechtedl 19761 ), which describes the LF with a power-law 
faint end and an exponential cutoff at the bright end. It 
is characterized by three parameters: M* , the characteris- 
tic bright-end magnitude, a, the faint-end slope, and $*, a 
normalization parameter: 

$(M) = (0.4;nlO)<I>*10 ' 4(M *- M)(1+a) e xp(-10 ' 4(M ^ M) )(9) 

As our method requires us to model the contribu- 
tions from LFs of all galaxy populations at all redshifts to 
the observed ga l axy d istribution simultaneously, we follow 
iGiallongp etall (|2005l ) in using an additional parameter S 
to model evolution in the parameter M*\ 

*(M,2) = (0.4 In 10) 10 0-4(M*-4io 910 [(l+,)/(l+ Z o)]-M)(l+ a ) 
Xexp t_lQ0A{M*-6log x0 [0-+z)/(l+z )]-M) 

This equation is derived from Eq. [9] simply by substi- 
tuting M* at with a redshift-dependent M*(z) = Mq — 
Slogio[(l + z) / (1 + zo)]. Here, Mq is the characteristic bright- 
end magnitude evaluated at a fiducial redshift zo. To mini- 
mize covariances between Mq and 5, zo should be chosen as 
the redshift at which M* is best-constrained by the sample. 
In our case, we have experimentally found that, for both 
early and late-type galaxies, the Mq-8 covariance is mini- 
mized for zo ~ 0.95. 

For modeling the LF as a function of SED type and red- 
shift, we assume that each SED template type is described 
by a Schechter function with a constant a, M* , and S. As 
photometric data cannot resolve density fluctuation on ar- 
bitrarily small scales, we assume that the LF for a given 
SED type is constant in each of ten redshift bins that we di- 
vide the redshift range 0.02 ^ z < 6 into. Potential density 
fluctuations on scales similar to or larger than these bins 
can therefore by reproduced by the algorithm. The redshift 
bins, along with the approximate limiting absolute magni- 
tudes for early and late type galaxies, are shown in Table [1] 
The divisions between redshift bins are motivated by the de- 
sire to apportion approximately equal numbers of observed 
galaxies to each bin. Also, we have attempted to reflect the 
robustness of photometric redshift constraints in the choice 
of binning. Although this model represents a compromise, 
the number of free parameters is still considerable. 

Finally, to provide a visual verification of whether the 
Schechter function is indeed an adequate description of the 
overall shape of the LF, we provide a binned representation 
of the LF. Here the free parameters are simply the averaged 
values of the LF in a number of discrete bins: 

${M) = Y;®iH{M,i), (11) 

i 

where H(M, i) = 1 if Mi < M < M i+1 and H(M, i) = 
otherwise. If evolution of the parameters M* and $* is sus- 
pected across the redshift range for which this representation 
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Table 1. Rcdshift bins 



bin# 


z 


yearly 
lira 


lim 


1 


0.02 - 0.1 


-11.5 


-11.5 


2 


0.1 - 0.25 


-14 


-14 


3 


0.25 - 0.4 


-15.75 


-15.5 


4 


0.4 - 0.6 


-17 


-16.5 


5 


0.6 - 0.8 


-18.5 


-17.75 


6 


0.8 - 1.2 


-20 


-19 


7 


1.2 - 1.5 


-21.5 


-20 


8 


1.5 - 2.5 


-22.75 


-20.75 


9 


2.5 - 3.5 


-23.75 


-21.75 


10 


3.5 - 6.0 


-24.5 


-22.75 



Broad rcdshift bins within which we analyze the LF. M, arly and 

J l%ra 

M|?*f give approximate effective magnitude limits around which 
the completeness of our sample drops below ~ 10% for average 
early- and late-type SEDs. 

is used, the interval limits Mi can of course be taken to be 
relative to the evolving M* at a given redshift, and an ad- 
ditional redshift-dependent normalization constant can be 
introduced to take out any previously known redshift evolu- 
tion of $* . This is the approach that we will be following in 

333 



3.4 Error Calculation 

The PML returns a likelihood parameter S = — 2 In C for 
every set of trial LF parameters. Therefore, in principle, 
uncertainty contour s on any pa rameter can be determined 
with the method of lAvnil (|l976h . i.e., the uncertainties are 
given by the set of parameter values of a constant likelihood 

= £-max X e ^ , (12) 

or, in terms of the likelihood parameter S in Eq. , 

S = S mtn + A x 2 - (13) 

Here, L ma x is the maximum likelihood, and A\ 2 is the 
\ 2 value for the desired level of confidence and the appro- 
priate number of degrees of freedom. The latter is usually 
the number of parameters whose uncertainties we are ex- 
amining, i.e., two if we are considering the uncertainties in 
the M* — $* plane. Following the example of Avni, we call 
these parameters the "interesting parameters", while refer- 
ring to all other parameters that the likelihood may or may 
not depend on, but whose values do not interest us at the 
moment, as the "uninteresting parameters". Eq. 1 13 1 requires 
that, for any choice of interesting parameters whose likeli- 
hood we wish to probe in order to find their uncertainty 
contours, the likelihood be maximized with respect to all 
uninteresting parameters. 

There are some parameters in our model for the LF and 
its evolution that we include neither in the best-fit search 
nor in the uncertainty calculation; we refer to these as back- 
ground parameters. In our case, we treat the LF parameters 
of the very high and very low redshift LFs, which are ill- 
constrained by our data and sensitive to systematic errors, 
as such background parameters. We set these to realistic val- 
ues motivated by extrapolating the a, M* , and $* values 
from the constrained redshift bins in order to represent the 



approximate contribution of the high- and low-redshift uni- 
verse to the observed galaxy counts. We do evaluate their 
impact on our fits with a Monte Carlo approach by adding 
large Gaussian perturbations to their values and repeating 
the best-fit search. This approach is not meant to repre- 
sent a realistic propagation of the uncertainties of these 
high-redshift parameters, because neither are the sizes of 
the propagated perturbations scientifically motivated, nor 
are these parameters necessarily statistically independent. 
Rather, this procedure demonstrates the covariances that 
exist between these background parameters and the param- 
eters constrained by our analysis, and serves to illustrate 
whether a given recovered LF parameter is robust against 
such implicit or explicit assumptions about the LF in red- 
shift regions where it cannot be constrained directly. 



3.5 SED Fitting 

The Photometric Maximum Likelihood method for the cal- 
culation of the LF, just like photometric redshift algorithms, 
requires a set of spectral energy distribution (SED) tem- 
plates for galaxies in order to map the assumed LFs into 
a galaxy distribution in photometric space. The templates 
we have chosen are based on th e empirical templates by 
IColeman. Wu fc Weedmanl (|l980l), extrapolated into the UV 
using Bruzual- Chariot models lIBruzual fc Chariot! |l993h . 
These base templates correspond to Ellipticals (E), Sbc 
galaxies, Scd galaxies, and very late-type irreg ulars (Im). In 
additi on, we use a starburst template (SB1) bv lKinnev et al.l 
( 1996). We create linear combinations of these templates and 
furth er apply a Calzetti dust extinction law l|Calzetti et al.l 
|2000T ) with a range of extinctions to them, arriving at a base 
template set comprising over one hundred SEDs. We then 
optimize this template set with regard to its ability to repro- 
duce the redshifts of the subsample of galaxies with available 
spectroscopic redshifts. The optimization is carried out by 
determining for each possible template whether its inclu- 
sion in the final template set would improve the quality of 
the best-fit photometric redshift solution with regard to the 
rms photometric redshift residual and the fraction of highly 
discrepant outliers. Templates are accordingly added to and 
removed from the set iteratively until no further improve- 
ment is possible. Table [5] describes how the 12 final tem- 
plates are generated from the four CWW templates. One 
template, "90% Im + 10% SB1 mod.", has been generated 
from a late-type one by strongly suppressing flux below 4000 

A. 

In the course of our analysis, we group several of these 
templates into two broad classes, early and late types. The 
possibility of a finer distinction exists, up to assigning in- 
dividual LFs to each SED. We have experimented with up 
to five distinct SED classes, but have found the recovered 
LF shapes to fall into two broad categories, one covering the 
two earliest-type templates, and one covering all later-type 
templates. Therefore, we will base our discussion of the LF 
results on just the two aforementioned types. In the present 
analysis, each individual template is weighted equally, i.e., 
every template belonging to the same class (early or late) is 
assumed to be described by the same LF, including the nor- 
malization. In subsequent work, we will explore the option of 
weighting templates differently in order to reproduce a more 
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Table 2. Grouping of SED templates into SED classes 



Template 



Class 



E 






early type 


80% 


E + 


20% Sbc 


early type 


Sbc 






early type 


Scd 






late type 


80% 


Scd 


+ 20% Im 


late type 


80% 


Scd 


+ 20% Im, E(B-V)=0.05 


late type 


40% 


Scd 


+ 60% Im 


late type 


30% 


Scd 


+ 70% Im 


late type 


10% 


Scd 


+ 90% Im 


late type 


Im 






late type 


90% 


Im - 


h 10% SB1 


late type 


90% 


Im - 


- 10% SB1 mod. 


late type 



The first column gives the linear combinations of the five base 
templates E, Sbc, Scd, Im, and SB1, that we use in our analysis. 
One template is slightly dust-reddened. The second column shows 
how templates are assigned to the two SED classes, early and late 
type- 
realistic representation of the distribution of the galaxy pop- 
ulation in photometric space. 

Although we describe the 12 templates using just two 
LFs, the large number of templates nonetheless has an im- 
portant effect: it broadens the footprint of a given LF in 
photometric space, allowing it to account for galaxies that 
may be incompatible with having been drawn from just any 
one individual SED. Table [2] shows how the templates are 
distributed into early and late types. 

An obvious caveat in this approach is that the empirical 
templates are appropriate for bright low-redshift galaxies, 
but may not be representative of the high-redshift universe, 
nor of very faint dwarf galaxies. Furthermore, the empirical 
optimization of our template set is naturally biased towards 
templates that represent the type of galaxies in the spectro- 
scopic subsample, which was not chosen systematically. For 
these reasons, we take special care to verify the robustness of 
our conclusions against assumptions about the high-redshift 
LF (and thus, by extension, about the SEDs that we asso- 
ciate these high-redshift LFs with). 

In subsequent studies of the full MUSYC sample, we 
plan to address these restrictions by allowing for a) evo- 
lution of the template set with redshift and b) superposi- 
tions of templates to overcome the discreteness inherent in 
the current template set. Both approaches can be integrated 
naturally into the PML, as long as an unambiguous assign- 
ment can be made between a certain SED template or tem- 
plate combination and the LF that is supposed to describe 
this type of galaxy. 



3.6 Verification of the PML with a Mock Catalog 

3. 6. 1 Construction of the Mock Catalog 

We have carried out extensive tests of our algorithm to 
verify its ability to recover the LF and other physical pa- 
rameters of the galaxies in our sample. Most of these tests 
are based on applying the method to a mock galaxy cat- 
alog. We use a mock catalogue extracted from a ACDM 
numerical simulation populated with GALFORM (version 
corresponding to Baugh et al., 2005) semi-analytic galax- 



ies. The simulation consists of a total of 10 dark-matter 
particles in a cubical box of 1000h _1 Mpc a side, followed 
from an initial redshift z — 30. The background cosmol- 
ogy corresponds to a model with matter density parameter 
£l m = 0.25, vacuum density parameter £7a = 0.75, a Hubble 
constant, H = 7il00kms _1 Mpc , with h = 0.7, and a pri- 
mordial power spectrum slope, n s = 0.97. The present day 
amplitude of fluctuations in spheres of 8h _1 Mpc is set to 
erg = 0.8. This particular cosmology is in line with recent 
cosmic microwave background anisotropy and large scale 
structure measurements (WMAP team, Spergel et al. 2007, 
Sanchez et al., 2006). 

The catalog provides, among other parameters, red- 
shifts, r-band luminosities, and bulge luminosity fractions 
(i.e., the fraction of the total luminosity that is associated 
with a bulge component of the galaxies). The LF in the 
simulation output is independent of redshift. The simula- 
tion output is volume- and luminosity-limited and complete 
to M r = — 16 and z=3. 

We have chosen a spectral energy distribution for each 
mock galaxy by superimposing up to four different SED tem- 
plates. The templat es used for this purpose are the fou r em- 
pirical templates bv lColeman. Wu fc Weedmanl (1980), two 
starburst temp lates by Kinney et al.l (1 19961 ) and the SDSS 
QSO template (jVanden Berk et al H200C The templates are 
picked randomly, with a probability distribution that varies 
depending on the bulge fraction of the input galaxy. There- 
fore, the range of possible SEDs for the mock galaxies far ex- 
ceeds those represented by the 12 SED templates used in our 
analysis, and the mock SED generally does not correspond 
exactly to any one of the templates used for the LF recovery. 
We then obtain expectation values for the galaxy photome- 
try in the UBVRIz J7f-bands by applying the transmission 
functions (filter curves, atmospheric transmission, and quan- 
tum efficiency) of the MUSYC survey, while simultaneously 
preserving the flux in the r-band specified in the output of 
the cosmological simulation. Uncertainties are taken from 
the MUSYC galaxy that best matches the expected fluxes 
under a least-squares comparison. 



3.6.2 Photometric Redshift Recovery 

The PML, just like any photometric redshift algorithm, is 
based on the assumption that the full range of physical char- 
acteristics of a galaxy (redshift, spectral type, absolute mag- 
nitude) can be recovered from the photometric properties of 
a galaxy. A simple way of demonstrating that the data are 
of sufficient quality and that our code is capable of recover- 
ing this information is to extract the information that would 
normally be provided by a photometric redshift analysis, i.e., 
constraining z and SED through the photometric vector. 
This is easily accomplished by finding, for each galaxy i, the 
model parameters zo, SEDo, and Mo for which the likeli- 
hood term p(fi | zo, SEDo, Mo) is maximized. This is equiv- 
alent to the best-fi t solution yielded by photo-z algo rithms 
such as HYPERZ jBolzonella. Miralles fc Pelldll2000h . 

Fig. [2] shows the recovered photometric redshifts ver- 
sus the true redshifts in the input catalogue for 23256 mock 
galaxies with mu < 25, split by the recovered best-fit SED 
template type (early or late). The figure shows a generally 
satisfying agreement between the input redshift and the re- 
covered photometric redshift, especially for early- type ob- 
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Figure 2. Best-fit photometric redshifts recovered by our LF 
algorithm from the mock galaxy catalogue. Upper panel: Galaxies 
best fitted with early-type templates. Bottom panel: Galaxies best 
fitted with late-type templates. 



jects. The photometric redshifts for late-type galaxies are 
less-well constrained, presumably due to the lack of charac- 
teristic structures such as the 4000 A Balmer break in early- 
type galaxy spectra. There are preferred ranges for the pho- 
tometric redshift, clearly visible as horizontal bands in the 
diagram, presumably indicative of redshifts at which charac- 
teristic breaks in the SED fall between neighbouring filters. 
It is worth pointing out that this is how purely Gaussian er- 
rors in the photometry propagate, i.e., assuming a Gaussian 
distribution for photometric redshift errors is, in general, a 
fallacy. 

The standard deviation of the redshift residual (z p hoto — 
ztrue) I (1 + ztrue) for all galaxies down to m — 25 is a = 0.40 
(including catastrophic outliers); when rejecting outliers of 
3a or more, the residual standard deviation is a — 0.059 
(with 13% of the objects having been rejected as out- 
liers). The accuracy is slightly better for the early types 
(a = 0.048) than for the late types (a = 0.061). There is 
a slight systematic offset of Az/(1 + z) = 0.005 (0.002 for 
early types and 0.005 for late types). 



3.6.3 Luminosity Function Recovery 

We will now demonstrate that our algorithm is able to re- 
cover the LF from the mock galaxy catalog by carrying out 
three different analyses of the mock catalogue. As the mock 
catalogue is populated using a semianalytic galaxy forma- 
tion model, and not a model Schechter function, we need 
to recover parameters for the "true" LF from the catalogue 
first. We do this by constraining mock galaxies to their true 
redshift and the best-fit SED and absolute magnitude at 
that redshift. In the second analysis, we will use the best-fit 
photometric redshift, SED, and absolute magnitude, equiva- 
lent to the conventional approach of substituting photomet- 



ric for spectroscopic redshifts. The third analysis will consist 
of applying the full PML formalism to the mock catalogue. 

The procedure for recovering the LF will be identical in 
all cases. However, in order to constrain a galaxy to a given 
redshift, SED type, and absolute magnitude, we modify Eq. 
[3 to read 

\i = / dM j dz \ $( Z — Ztrue) 

SED ** ** 

x6{SED - SED t rue)5(M - Mtme) x <P(Mo;SED;z | P), (14) 

where ztrue, SEDtrue, and Mtrue are the values to which a 
given galaxy is constrained; th i s app roach is similar to the 
original one by iMarshall et aL (1983). In the calculation of 
the second term in Eq. [8] we calculate A(/) as previously 
(an integral over all of observable space), but apply a hard 
cutoff at the survey magnitude limit without allowing for 
photometric uncertainties. This procedure is then equiva- 
lent to running the Marshall method on a parameter space 
comprising z, M, and SED, rather than observed fluxes. 

For each of the three analyses, we then follow the pro- 
cedure that we will subsequently apply to the real data: We 
model the LF as a function of two SED types (early and 
late) in ten redshift bins. Over the range 0.1 ^ z < 1.2 (red- 
shift bins 2 to 6), we describe it in terms of Eq. 1101 with a 
constant faint-end slope a and a two-parameter model for 
M* and its evolution with redshift for each SED type. The 
normalization parameters <&* are independent in all bins; 
<E>* are included in the optimization in redshift bins 7 and 8. 
All other LF parameters, specifically in the lowest and high- 
est redshift bins, are held constant during the optimization, 
but periodically readjusted to match extrapolations from the 
constrained bins. 

An added complication is that the early-type LF in the 
input mock catalogue is not a single Schechter function, but 
exhibits a possible bimodality and a sharp cutoff at M r = 
— 16, neither of which is represented by the LF models that 
we are fitting. However, this complication will affect all three 
analysis methods. 

Fig. [3] shows the 68% and 95% error contours in the a- 
Mq plane for all three analyses: the calculation based on the 
true redshifts in dotted lines, the calculation using best-fit 
photo-zs in solid thin lines, and the PML solution in thick 
lines. Both for early and late type populations, the conven- 
tional approach, using best-fit photo-zs, fails to recover the 
correct LF parameters in the M* — a plane. Although the 
fits would appear consistent if projected only onto the M* 
axis, a comparison involving only M* must be carried out 
at the same a, in which case the discrepancy between the 
photo- z-based solution and the true LF is even more appar- 
ent. 

The late-type LF parameters recovered by the PML, 
on the other hand, are consistent with the solution obtained 
from the true redshifts. For the early-type population, there 
is a discrepancy in a, so that the 2a contours barely over- 
lap; however, for an a ~ —1.25, which is consistent with 
both the PML analysis and the true redshift calculation, 
the recovered M* are also perfectly consistent, unlike the 
conventional, photo-z-based solution, which is significantly 
brighter. 

The PML solution can be brought into even better 
agreement for different choices of values for the constant 
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background parameters at very low and high redshift. We 
have carried out the procedure described in ^3.41 to deter- 
mine how uncertainties in the LF priors in the very high- 
and low-redshift bins propagate into the LF fit. For each 
trial calculation, we apply Gaussian perturbations of order 
1 mag to the M* background parameters (i.e., those not in- 
cluded in the optimization) , and of order unity to the quan- 
tity logio<&* (i.e., M* typically scatter randomly within 1 
mag, and $* within a factor of 10 from the default values, 
which are very generous assumptions). We then calculate 
the best-fit solution for the constrained parameters. 

The best-fit results of these Monte Carlo trials are 
shown in Fig.[3]with dots (for all fits) and crosses (for those 
fits that yield an equivalent or better likelihood after opti- 
mization than our default scenario). Each symbol represents 
one Monte Carlo trial with a different LF prior in the low- 
and high-redshift bins, after the maximum likelihood opti- 
mization of the free parameters has been carried out. 

A large fraction of the Monte Carlo trials yield a better 
total likelihood than our default fit. Since we do not apply 
any prior to judge which of the various background LF pa- 
rameters that are probed by the Monte Carlo algorithm are 
realistic, and because, for reasons discussed in greater detail 
in £1-4.31 we do not consider our data set to be ideally suited 
for constraining the high-redshift LF, we do not actually 
prefer any of the Monte Carlo trials with better likelihood 
to our default fit; however, regardless of whether the back- 
ground LFs assumed for the Monte Carlo trial is realistic or 
not, runs that yield a best-fit likelihood that is significantly 
worse than our default scenario can be ruled out; therefore, 
only the crosses in Fig. [3] should be considered. 

We find that, of the Monte Carlo trials that result in an 
equivalent or better likelihood than the default fit, almost 
all yield LF parameters that are in even better agreement 
with the true LF; in fact, for some background priors, we 
recover LF parameters that are in almost perfect agreement 
with the input LF (we have not introduced criteria, however, 
to decide whether the background LF priors corresponding 
to these trials are physically plausible or not, which is why 
we do not elect any of these as the real best fit). Conversely, 
of those LF parameters that are in worse agreement with 
the true LF than the default fit, almost all are marked by 
having a substantially worse likelihood and are thus ruled 
out. This means that the PML still has some constraining 
power in redshift regions that we have chosen not to include 
in our optimization in our conservative approach. 

However, the Monte Carlo approach does not open up 
additional degrees of freedom to the photo-z-based algo- 
rithm, because (at least if the photo-zs are derived from 
frequentist best-fit estimates), each galaxy is uniquely asso- 
ciated with a given redshift bin, independently of the LF 
assumed for it. Therefore, the photo-z-based LF remains 
clearly inconsistent with the true LF. 

The advantage of the PML method over the use of pho- 
tometric redshifts is even more obvious if we reduce the num- 
ber of filter bands. Using only the BVRIz filters, we recover 
the results shown in the inset in Fig. [3] Here, the photo-z- 
based results are drastically different from the "true" LFs; 
the late-type LF is found to be a ~ —1.3, compared to 
the actual a = —1.52, and the recovered early-type LF is 
also significantly brighter and shallower than the input LF. 
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Figure 3. Uncertainty contours at the 68% and 95% levels for 
the LF parameters recovered from the mock catalogue using true 
redshifts (dotted contours), best-fit photo-zs (thin solid contours) 
and the PML algorithm (thick solid contours). The main panel 
shows the results from on UBVRIzJK filter photometry, the in- 
set only BVRIz. Dots and crosses demonstrate how fluctuations 
in the high- and low-redshift LF parameters not included in the 
optimization may propagate into the recovered LF parameters, 
as described in the text (crosses symbolize a likelihood equal or 
better than the default scenario). The photo-z-based LFs are in- 
consistent with the true LF, while the PML solution overlaps with 
the true LF at at least the 2a level. The distribution of crosses 
shows that different LF priors in the lowest- and highest-rcdshift 
bins might bring the LF solutions in almost perfect agreement. 
The performance differences between photo-z-based calculations 
and the PML are particularly striking for the smaller filter set. 



In contrast, the uncertainty contours of the PML solution 
again overlap with those of the "true" LF. 



4 RESULTS AND DISCUSSION 
4.1 Magnitude Limits 

How deep does MUSYC probe the galaxy LF in our individ- 
ual redshift bins? To answer this question, we plot in Fig. 
[4] the derivative of the second term of the likelihood func- 
tion in Eq. [8] with respect to the value <E>(M) of the LF at a 
given absolute magnitude; this is, in essence, the volume over 
which a galaxy with a given SED and absolute magnitude in 
a given redshift bin would be visible. It has here been nor- 
malized to unity for the brightest galaxies. The figure shows 
the completeness for the earliest- and the latest-type SED 
in redshift bins 2 (0.1 ^ z < 0.25), 3 (0.25 ^ z < 0.4), 4 
(0.4 < z < 0.6), 6 (0.8 2 < 1.2), and 8 (1.5 < z < 2.5). 

For higher redshifts, the early- and late-type complete- 
ness curves differ significantly, indicating that late-type 
galaxies are more easily detectable than early-type ones. The 
curves do not cut off abruptly at a given absolute magnitude 
due to the fact that each redshift bin has a finite extent in 
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Figure 4. Effective r-band magnitude limits for five of our ten 
redshift bins. Different line styles indicate different redshift bins. 
For each redshift bin, we draw two completeness curves, corre- 
sponding to very late and very early types. Assignment of line 
styles to redshift ranges is given in the plot. 



Figure 5. Best-fit photometric redshifts compared to known 
spectroscopic redshifts in the MUSYC catalog. Upper panel: 
Galaxies best fitted with early-type templates. Bottom panel: 
Galaxies best fitted with late-type templates. 



redshift; galaxies near the low-redshift boundary of a bin are 
visible to fainter absolute magnitudes than those at the high- 
redshift boundary. In the lowest redshift bin under analysis 
here, our sensitivity extends to Mr » —14. 

4.2 Photometric Redshift and Luminosity 
Recovery 

We have already demonstrated that the PML is capable 
of recovering the redshift of galaxies in the mock catalog. 
In the application to real data, additional systematic er- 
rors that are not modeled by the mock catalog may come 
into play, such as systematic biases in the photometry that 
our zero-point corrections may not have sufficiently com- 
pensated for. The best-fit photometric redshifts recovered 
by our algorithm are a convenient way of testing for the 
presence of such serious systematic errors. For this compar- 
ison, we have drawn spectroscopic redshifts in the E-CDFS 
field from various public and proprietary sources, including 
the NASA Extragalact ic D atabase (NED) , publ ications by 
iRavikumar et all (j2007r i and lVanzella et all (2008), as well as 
redshifts determined within the MUSYC project, yielding a 
total of 2602 redshifts. 

Fig. [5] shows the performance of the recovered best- 
fit photometric redshifts versus spectroscopic redshifts. The 
plot resembles our previous comparison using the mock 
catalogues in Fig. [2] The performance is generally good; 
late-type objects suffer from more catastrophic failures. 
The standard deviation of the redshift residual (z p hoto ~ 
Ztrue)/{1 + ztme) is (T — 0.27 (including catastrophic out- 
liers); when rejecting outliers of 3er or more, the residual 
standard deviation is a — 0.062 (with 17% of the objects 
having been rejected as outliers). The accuracy is similar 



for early types (a = 0.061) and for late types (a — 0.063). 
For comparison, prior to applying our photometric correc- 
tions to the zero-points and uncertainties, these numbers are 
a = 0.32 and a = 0.082, respectively (a = 0.077 for early 
types and a — 0.083 for late types). 

Comparison between Fig.[2]and Fig.[5]shows that many 
of the structures (particularly the horizontal patterns show- 
ing preferred values for photometric redshifts) found in the 
real mock data are visible in the observed data as well. This 
provides important verification that the pattern of catas- 
trophic failures in the photometric redshift recovery can be 
reproduced on the basis of just the Gaussian errors in the 
input photometry and that, by extension, such color degen- 
eracies that normally present a grave problem for photo-z 
based LF determines can be accounted for naturally by the 
PML. 

However, there is an indication of photo-zs of early type 
galaxies being systematically underestimated at z > 1; this 
feature is only partially reproduced by the mock data, and 
thus not fully corrected for by the PML either. Therefore, we 
restrict our analysis to the redshift range 0.1 ^ z < 1.2, the 
last redshift bin covering the fairly wide range 0.8 ^ z < 1.2. 

The fact that our photometric data and SED template 
set allow us recover photometric redshifts with acceptable 
accuracy is encouraging, even though our algorithm does not 
use individual photometric redshifts per se. For the purpose 
of calculating luminosity functions, the ability to recover the 
absolute magnitude of a galaxy and account for fc-corrections 
and other cosmological effects is at least as important. Al- 
though, unlike the redshift, the absolute magnitude cannot 
be verified by alte rnate observations, the COMBO-17 survey 
llWolf et aLll2004l) offers a high-quality comparison catalogue 
that covers the same field and contains estimates of absolute 
magnitudes in the r-band. Crossmatching the COMBO-17 
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catalogue and the MUSYC catalogue yields 26234 galax- 
ies. After adding an offset of 0.055 mag to the COMBO- 
17 R-band magnitudes to con vert from the Vega to the 
AB system IjFrei fc Gunnlll994f ). the apparent magnitudes 
in both catalogues agree well apart from a systematic offset 
of rriMusYC — rnaoMBO-n ~ +0.045 (in the .R-band; total 
magnitudes for COMBO-17, and either corrected aperture 
or AUTO magnitudes for MUSYC, as described in £[2j) . 

To evaluate our ability to recover absolute magnitudes 
from the MUSYC photometric data, we now constrain the 
MUSYC galaxies to be at the photometric redshift listed in 
the COMBO-17 catalogue, and carry out an SED fit at that 
redshift to calculate the r-band rest frame absolute mag- 
nitude. For a given redshift and SED type, the PML cal- 
culates the rest-frame luminosity of a galaxy directly from 
the normalization factor that yields the best fit of the SED 
to the observed fluxes. We compare thes e absolute magni - 
tudes to the r-band estimates provided by I Wolf et al.l ( 2008) 
by calculating the quantity (Mmusyc — MooMBO-n) ~ 
{mMUSYC — mcoMBO-n), which accounts for the mea- 
sured differences in observer-frame R-band apparent magni- 
tudes to isolate only the differences in the recovered absolute 
magnitude. Offsets between the A B and Vega system s are 
negligible in the rest-frame r-band (|Frei fc Gund[l994r ) . We 
restrict this comparison to galaxies with R < 22, because 
fainter galaxies (R > 22.5) in the COMBO-17 catalogue ex- 
hibit a strong quantization in the best-fit photometric red- 
shifts that casts doubts on their reliability. We then find a 
mean of ~ —0.06 with a dispersion of 0.22 over the red- 
shift range 0.1 ^ z < 1.0, i.e., MUSYC galaxies are inferred 
to be slightly more luminous at the same redshift and for 
the same rest-frame R-band magnitude. This value is only 
very weakly dependent on redshift (with a slope of —0.05 
and a zero point of —0.04 at z=0). A small contribution 
to this discrepancy may come from differences in the cos- 
mological model in the calculation of absolute magnitudes 
(Q\ = 0.73 and f2 m = 0.27 in our study, Q\ = 0.7 and 
fi m = 0.3 for COMBO-17); this difference is ~ 0.04 mag at 
2 = 1 in the sense of MUSYC galaxies being inferred to be 
brighter. Therefore, if residual systematic discrepancies exist 
in our calculation of absolute rest-frame magnitudes in com- 
parison to COMBO-17, they are at the level of < 0.1 mag for 
the redshift range that we are analyzing. Possible reasons for 
this residual discrepancy include the fact that our rest-frame 
luminosities are constrained by all available filter bands, so 
that the subtraction of the {mMUSYC — mcouBO-vi) term 
in the R-band will not account for all discrepancies in the 
input photometry, as well as the different template set. 

4.3 The Cosmic LF 

We now apply the PML to the MUSYC data to calculate the 
field galaxy LF. For this purpose, we use the parametriza- 
tion given by Eq. 1101 i.e., for each of the two SED classes, 
we assume a characteristic bright-end magnitude Mq at 
a fiducial redshift z = 0.95 (experimentally found to be 
the redshift for which M* is constrained best in this sam- 
ple), a faint-end slope a, and an evolutionary parameter S 
to describe the variation of M* with redshift according to 
M* — Mq — S x lg[(l + z)/(l + z )]. The normalization pa- 
rameters $* are independent for each redshift bin and SED 
class. For the present work, we have chosen to calculate the 



LFs in the rest-frame r-band. This band is less sensitive 
to effects of current star formation than, for example, the 
B-band, but at the same time is well-constrained by the 
available photometry over the redshift range of particular 
interest to us. 

The purpose of this calculation is to recover an approx- 
imate representation of the LF, albeit with a low number of 
free parameters, which will serve as the starting point for a 
subsequent refinement with a more detailed ansatz. Further- 
more, we will use this calculation to constrain the faint end 
slopes a for the early and late SED types from this calcu- 
lation and adopt them for the rest of the analysis. While it 
would be possible to optimize with respect to the faint end 
slopes a in individual redshift bins, the flux limit prevents 
us from obtaining useful constraints on a in any but the 
lowest-redshift bins, up to z m 0.4 — 0.5. A follow-up study 
of a larger sample will explore this issue. 

Parameters that we include in the optimization are Mo, 
a, and S for both SED types, describing the shape of the LF 
over the range 0.1 ^ z < 1.2. Furthermore, we constrain $* 
over the range 0.1 ^ z < 2.5. All other LF parameters (i.e., 
a, M* in redshift bins 1 and 7 to 10, and $* in bins 1, 9, 
and 10) are held constant. The decision to treat the high- 
redshift LF parameters as constant background parameters 
is based on a number of criteria: 

• The rest-frame r-band is not constrained directly by the 
available filter complement at high redshift. 

• Only the tip of the LF is sampled, which may introduce 
strong systematic errors if the most luminous galaxies do not 
obey the same Schechter function implied for the bulk of the 
galaxy population. 

• The applicability of z=0 empirical templates to high- 
redshift galaxy populations is doubtful. 

• Photometric redshifts, using this catalogue and our 
template set, have proven much less reliable at z > 1.0, 
as is also shown by Fig. [5] and Fig. [5] 

At low redshift, on the other hand, constraints on the 
LF are weak due to the small number of galaxies and the 
large uncertainties affecting the absolute magnitudes of in- 
dividual galaxies. 

The initial values that we choose for t hese background 
param eters are motivated by the results of lGiallongo et ail 
(2005). Their work uses two different criteria to separate 
early and late type galaxies — colors and specific star forma- 
tion rate — , and examines two ways of defining the threshold 
values — empirically, from the observed bimodality, or from 
a theoretical evolutionary track. This yields four different es- 
timates for the LF and its evolution for a given galaxy type. 
To obtain an estimate in a given redshift bin, we evaluate 
their fit at the central redshift of the bin, and average the 
four extrapolated values for M* and $*. To convert their B- 
band LFs into our r-band rest-frame absolute magnitudes, 
we adjust M* by B — r = 1.35 for early types and B — r = 0.9 
for late types. These values corresp ond to typical B— r colors 
for early and late populations from I Wolf et al.l (|2003T > at low 
redshift. We caution that the distinction made by Giallongo 
et al. between early and late types is different from ours; in 
particular, while we use redshift-zero SED templates to fit 
early and late type LFs, their threshold values for parame- 
ters such as color and specific star formation rate vary with 
redshift. We therefore expect that "early" types would be 
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more prevalent at high redshift according to their definition 
than to ours. 

We repeat the best-fit search several times; after each 
iteration, we adjust the background parameters by averaging 
their previous values with linear extrapolations of the best- 
fit M* and values over lg(l + z). Above the last redshift 
bin directly constrained by our fit (above z—2.5), we cap the 
extrapolated values for these parameters at those found in 
the closest redshift bin where they are constrained directly, 
so as to avoid implying unrealistically bright or abundant 
galaxy populations in the high-redshift universe. 

We solve for the maximum-likelihood solution in the 
remaining 16-parameter space (2 SED types x (3 shape 
parameters + 7 normalization parameters)) using an iter- 
ative algorithm that, starting from the previous solution, 
first searches for the maximum likelihood point along each 
parameter axis and thus deriving a shift Apt in each parame- 
ter pi. Subsequently, it applies a shift xAp to all parameters 
simultaneously and calculates the likelihood as a function 
of the scaling factor x in the same manner. If the resulting 
adjustment in any given parameter is much smaller than the 
spacing between the trial points that are initially evaluated 
along the parameter axis, this spacing is slightly reduced. 
The algorithm stops when the trial point spacings along all 
parameter axes are much smaller than the typical statistical 
uncertainties. 

The best-fit values for the 16 parameters included in 
the optimization are given in Table [3] along with the val- 
ues adopted for the constant "background" parameters. The 
projection of the a-Mg plane is shown in Fig. HJ] along with 
the 68% and 95% error contours. Here, bold error contours 
indicate the results of the PML analysis, while thin lines 
represent the classical calculation using fixed photo-zs; the 
uncertainty regions are clearly disjunct. 

The qualitative behaviour of the fits is consistent wit h 
many prev i ous lit erature results (i.e.. lBromlev et alJ (|l99Sf) ; 
IWolf et ail (|2003T )) regarding the late- and early-type LFs, 
whose most important characteristics are a rather shallow 
slope for early-type galaxies and a significantly steeper 
slope for late-type galaxies. However, the steepness of the 
slope of our late-type LF is remarkable indeed; we find 
— 1.62lg'g2, compared to — 0. 62^g'o5 for the early types. 
Literature values for the total LF tend to prefer values 
of a = —1.2 to —1.3. However, it is also k nown that the 
late- ty pe LF is steeper than this average, and lBlanton et ail 
(2005) find a — —1.3 in the SDSS, but speculate (based on 
the fact that the SDSS is likely to miss many low-surface 
brightness objects) that it may be as steep as —1.5 or 
even steeper. Such steep slopes at the faint end can be 
easily reconciled with older values of a ~ —1.2 by noting 
that the bright end of the LF is dominated by early-type 
galaxies, and that the late-type LF with its steep slope 
will begin to dominate the total LF at several magnitudes 
fainter than M* , i.e., only surveys combining very deep 
samples with sufficient statistics to constrain the extreme 
faint end will be able to identify this as an upturn. It is 
conceivable that this dichotomy may also be identified as a 
"dip" in the overall LF, especially in denser environments, 
where the steep late-type LF is likely to have been depleted 
by environmental processes and its normalization $* 
lowered. Indeed, claims of such a "dip" feature in the 
overall LF have been raised in the past based on surveys 
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Coma cluster (|Biviano et al.l ll996). 

The early and late type populations show a very differ- 
ent evolutionary behaviour with redshift. For the early- type 
sample, we find 8 ~ 0, i.e., no evidence for luminosity evolu- 
tion up to z = 1.2. In contrast, the late-type sample is fitted 
by an evolving LF with 8 ~ 2.5, i.e., M* brightens by about 
0.75 mag from z = to z = 1. However, the uncertainties 
on the 8 parameters are very large, and even the late-type 
sample is formally consistent with very weak luminosity evo- 
lution. 

Our two-parameter fit to M*(z) can be extrapolated 
to z = to recover the local LF parameters; however, as 
M* is primarily constrained by galaxies at z !*s 0.95, uncer- 
tainties in 8 propagate strongly into the extrapolated value. 
Formally, we find M*(z = 0) = -21.854+° 2 2 ° Q ° for the early 
type and M*(z = 0) = -21.632t°^° for the late type pop- 
ulation. However, better constraints on the LF at z=0 can 
be obtained by fitting the LF parameters directly in indi- 
vidual low-redshift bins, and we refer the reader to the next 
section, ij4.4l where we will discuss this question in greater 
detail. 

As we have discussed in H3.41 we presume LF parame- 
ters at very low and very high redshift to be constant instead 
of constraining them directly from the data, as they are very 
susceptible to systematic errors. However, as LF parameters 
in different redshift bins are usually covariate, the exact LF 
solution recovered here depends on our assumptions for the 
values of these background parameters. We have applied the 
Monte Carlo approach described in £13.41 to gauge the sensi- 
tivity of the best-fit parameters to these background values. 
For each trial calculation, we apply Gaussian perturbations 
of order 1 mag to the M* background parameters, and of or- 
der unity to the quantity logio&* in these background bins 
(i.e., M* typically scatter randomly within 1 mag, and $* 
within a factor of 10 from the default values) . We then calcu- 
late the best-fit solution for the constrained parameters. The 
results are shown in Fig. |5]in the M^-a plane. Dots represent 
Monte Carlo runs for which the final, best-fit likelihood is 
substantially worse than for our default values; this means 
that the corresponding configuration of background param- 
eters is ruled out by the data, and the run can be discarded. 
Crosses represent those runs for which we achieve a similar 
or better maximum likelihood than for our default back- 
ground choices; these represent background configurations 
that cannot be ruled out by the data, although we do not 
judge whether the configuration of background LFs assumed 
in each Monte Carlo trial is realistic, and therefore do not 
prefer any of these solutions to our default scenario. Accord- 
ing to Fig. [6] variations in the background parameters may 
then introduce a scatter that exceeds the 2a uncertainty con- 
tours. However, the qualitative picture remains unchanged; 
in particular, none of the acceptable Monte Carlo trials yield 
a > —1.5 for the late-type population. Furthermore, only 
very few outliers in the Monte Carlo analysis yield LF pa- 
rameters consistent with the photo-z-based calculation for 
the late-type population, and none of them are consistent 
with the photo-z-based early-type LF. We furthermore point 
out that the amount by which the background parameters 
are allowed to scatter from the default assumption is very 
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generous. We interpret this as indicating that, to change 
the faint-end-slope a of the late-type population to values 
shallower than —1.5, $* at high redshift would have to be 
permitted to deviate by more than a factor of 10 upward 
from the values found at lower redshift, or M* by more than 
a magnitude towards the bright end, i.e., the high-redshift 
galaxy population would have to be much brighter and more 
abundant than at z < 1.2, which is a scenario that we con- 
sider unrealistic. 

In order to test the effect of excluding detections with 
questionable star/galaxy classifications, we have carried out 
an alternate fit including only objects with S/G ^ 0.8 (as 
opposed to our default of 0.97). The resulting best-fit val- 
ues differ slightly from those obtained above; most notably, 
Mq for the early types is fainter at Mg — —21.62, and the 
faint-end slope is slightly shallower at a — —0.43. However, 
there is virtually no change in the late-type LF, where the 
best-fit values are Mq — —21.66 and a — —1.61. As a and 
Mq are correlated, it is not immediately clear whether the 
shift in the early type best-fit parameters is due to a change 
at the bright or the faint end; however, as there is no obvi- 
ous tendency for objects with S/G > 0.8 to be associated 
with unusually high luminosities, we suspect that the change 
is due to the removal of compact sources with photometric 
properties characteristic of low-redshift early-type objects. 
We return to this question once more in i]4.7l where we anal- 
yse the shape of the LF in a binned representation. 

Just as doubts regarding the appropriateness of our em- 
pirical SED templates prevent us from constraining the LF 
in the high-redshift universe, we may also wonder whether 
a mismatch between our late-type templates and the pho- 
tometric properties of dwarf galaxies may have steepened 
the observed faint end slope artificially. Dwarf galaxies, typ- 
ically being of low metallicity, may have significantly dif- 
ferent spectral properties than the brighter objects that 
the CWW templates were originally drawn from. However, 
we would expect template insufficiency to lower the recov- 
ered LF value associated with the galaxy populations that 
are not represented well by the template, while the LF of 
other populations with similar photometric properties would 
be boosted instead. Therefore, a greater concern regard- 
ing the steep faint end lies with the question whether the 
high-redshift LF is represented well by our templates. The 
Monte Carlo procedure described above, however, demon- 
strates that the steep faint end slope is not very sensitive to 
the LF values assumed for the high-redshift population and 
cannot be pushed to values shallower than a ~ —1.5 even 
for the rather large background fluctuations assumed there. 

4.4 The LF in individual redshift bins 



Our best-fit solution described in £14.31 has been based on 
a relatively low number of free parameters. One of the dis- 
advantages incurred by this is that the constraints on M* 
are weak, and the extrapolation to redshift zero via the 
two parameters Mq and 8 is potentially affected by large 
errors. Furthermore, our assumptions (such as the simple 
two-parametric model for M* as a function of redshift) can- 
not be verified explicitly within the confines of this simple 
approach. 

To explicitly test these assumptions and obtain more 
accurate constraints on the LF at low redshift, we now re- 
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Figure 6. Uncertainty contours for the best-fit LF parameters at 
68% and 95% confidence levels in a-M* parameter space. Bold 
contour lines show the PML solution, while thin contour lines 
indicate a solution based on a conventional algorithm using best- 
fit photo-zs. The results reflect a dichotomy between the early- 
type LF with a shallow faint-end slope and the late-type LF with 
an extremely steep slope of a m —1.6. Dots and crosses again 
demonstrate how fluctuations in the high- and low-redshift LF 
parameters propagate into the recovered LF parameters. Only 
trial runs symbolized with crosses yield an acceptable likelihood. 



fine our solution by maximizing the likelihood with respect 
to the parameters M* and $* in individual redshift bins. 
We iterate both parameters for each SED type and for each 
redshift bin. Other than that, the same parameters are in- 
cluded in this optimization as in the previous section: $* in 
redshift bins 7 and 8, and the faint-end slope a in bins 2 to 
6 (assumed independent of redshift). 

The results of this analysis — in an A/*-$* projection 
— are shown in Fig. [7] Bold contours represent early types, 
and thin contours late types. Error contours are calculated 
for two degrees of freedom and 68% confidence. Table[4]gives 
the best-fit values in individual redshift bins. 

The early-type galaxies exhibit relatively little scatter 
in M* — all error contours lie in the range —22.2 < M* < 
—21.5. Given the relatively large uncertainties, there is no 
obvious trend of M* with redshift, confirming the conclusion 
that S m 0. The scatter in $* is larger than the uncertainties; 
in the seven redshift bins included in our analysis, $* scat- 
ters within a factor o f ~ 3. The sc a tter m ay be affected by 
large scale structure. iPhleps et al] l|2007h report an under- 
density at 0.25 < z < 0.4, which coincides with our very low 
data point in redshift bin 3. On the other hand, Ed ami et all 
( 2005), based on spectroscopic data, note several structures, 
including massive galaxy groups, at z « 0.7, which falls into 
our redshift bin 5, which is significantly elevated in compar- 
ison to the neighbouring bins. Our LF results are therefore 
qualitatively consistent with large scale structure features 
in the E-CDFS known from spectroscopic studies. The scat- 
ter introduced by these features is too large to permit any 
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Table 3. LF Results at 0.1 < z < 1.2 



bin 


z 


a 


M* 


S 


lg $* 




early types 


-0 62+ 04 


-21.885_ 075 


107+ - 780 

u,lu ' -0.720 




1 


0.02 - 0.10 


-0.639 


-21.918 




-2.502 


2 


0.10 - 0.25 








9 *?91 +0-043 
— 037 


3 


0.25 - 0.40 








9 79C+0.044 
z - ,zu -0.039 


4 


0.40 - 0.60 








9 KQ9+0.034 
^• UOZ -0.026 


5 


0.60 - 0.80 








9 4OO+0.025 
z -* oz -0.020 


6 


0.80 - 1.20 








9 7 c; fi +0.028 
Z -' OD -0.027 


7 


1.20 - 1.50 


-0.639 


-21.867 




-2 709+ 053 
z -' ua -0.056 


8 


1.50 - 2.50 


-0.639 


-21.849 




9 fiS^+ ' 077 


9 


2.50 - 3.50 


-0.639 


-21.828 




-2.850 


10 


3.50 - 6.00 


-0.639 


-21.802 




-2.953 




late types 


-1 62+ 02 
±DZ -0.02 


-22 345+ - 075 
zz.o'io_ Q 1QQ 


2 484+ 1 260 




1 


0.02 - 0.10 


-1.664 


-21.683 




-3.158 


2 


0.10 - 0.25 








-2 Q94+ - 025 


3 


0.25 - 0.40 








-3 U4+ 029 
"5- 14 *-0 .064 


4 


0.40 - 0.60 










5 


0.60 - 0.80 








-3 029+ 013 
o.uz»_ 064 


6 


0.80 - 1.20 








o -, 9 o+0.018 
o.izo_ 067 


7 


1.20 - 1.50 


-1.664 


-22.749 




-3.68»_ Q85 


8 


1.50 - 2.50 


-1.664 


-22.967 




-2 769+ - 008 


9 


2.50 - 3.50 


-1.664 


-22.967 




-3.252 


10 


3.50 - 6.00 


-1.664 


-22.967 




-3.280 



Results from the fit of SED-specific Schechter functions to the redshift range 0.1 ^ 
z < 1.2. In the first row of each section of the table, we give the shape parameters 
describing the best-fit LF in this range. Normalization parameters <&* are listed sep- 
arately for each redshift bin. Values without uncertainties are constant background 
values and not optimized. The parameter M* is evaluated at z = 0.95. 



reliable conclusion as to whether $* exhibits a trend with 
redshift. 

The uncertainty contours for the late-type LFs are af- 
fected by the strong covariance between M* and 3>* in the 
case of a steep faint end. The best-fit values form an approx- 
imate sequence, with the faintest, but highest-normalization 
fits being found at low redshift, while the high-redshift LFs 
are brighter by more than two magnitudes in M* , but also 
exhibit a lower $* (by a factor of 2 — 3); this is also consistent 
with our earlier finding of a positive evolution parameter S. 
Comparison of the corresponding Schechter profiles shows 
that the variation of the late type LF between redshift bins 
occurs primarily at the bright end, whereas the normaliza- 
tion of the LFs in this redshift range is almost constant with 
redshift at magnitudes M r > —20. However, redshift bin 3 
deviates significantly towards lower normalizations from the 
sequence along which the uncertainty contours for the other 
bins are aligned, which is again consistent with the claim of 
an underdensity in the E-CDFS in this redshift range. 

We have again determined the impact of varying the 
high- and low-redshift LF parameters held constant dur- 
ing these fits using a Monte Carlo procedure as described 
in t )3.4l The range of best-fit parameters is shown in Fig. 
[7] by dots for all trials, and crosses for trials yielding ac- 
ceptable likelihood. For the early-type LFs, shifts towards 
slightly brighter M* /lower $* are possible, but shifts in the 
opposite sense (fainter M*) are ruled out. For the late-type 
sample, the background uncertainties propagate primarily 



into variations in We remind the reader that the mag- 
nitude of these variations is dependent on the magnitude 
of the high- and low-z background uncertainties, which are 
assumed to be very generous here (a factor of 10 in $* , and 
one mag in M*), and that the Monte Carlo procedure can 
therefore only convey a qualitative impression of the robust- 
ness of the results and the sense in which they may deviate 
from the default solution. 

4.5 Comparison to other surveys 

Our results for the low-redshift LF appear to be in contrast 
to ot her determinations for the local field LF (|Blanton et alj 
2003); in particular, M* for the early- type population at 
low redshift appears strikingly bright, compared to a value 
of M* = -21.18 ± 0.01 at z=0.1 given by Blanton et al. 
(adjusted to Ho = 71). 

To illuminate the causes of this discrepancy, we present 
in Fig. [8] a comparison between the Blanton et al. LF at 
z=0.1 and our results in four low redshift bins, as given in 
Table|4] For the purpose of facilitating comparison, we have 
drawn the z — 0.1 Blanton et al. LF in all four panels, but 
emphasize that even the lowest-redshift bin of our study 
shown here lies above the redshift for which that LF was 
calculated. 

The histograms show the distribution of absolute mag- 
nitudes of galaxies with known spectroscopic redshifts in 
the MUSYC catalog. For bright galaxies, more than 25% 
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Table 4. LF Results at 0.1 < z < 1.2 
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Figure 7. Uncertainty contours at the 68% confidence level for 
fits to the early-type (bold contours) and late-type (light con- 
tours) LF in individual rcdshift bins over the range 0.1 z < 1.2. 
For the early-type LF, we see no evidence for an evolution of M* , 
while <!>* exhibits considerable scatter, which may be associated 
with known density inhomogeneities in the E-CDFS. The late- 
type LFs exhibit a strong covariance between M* and <E>* and 
populate a sequence ranging from faint Af*/high <E>* at low to 
bright M*/low $* at high redshift, attesting to evolution in at 
least one parameter. Dots and crosses illustrate how variations in 
the constant high- and low-redshift LF parameters may propagate 
into the recovered LFs. 



of the objects in the catalogue have spectroscopic redshifts; 
we have applied an approximate completeness correction as 
a function of the apparent magnitude ma and divided by 
the comoving volume of the redshift bin. Apart from the 
calculation of the rest-frame r-band absolute magnitudes, 
these are the only corrections applied to the data. In partic- 
ular, we have abstained from applying a volume correction 
in order to minimize the risk of introducing systematic er- 
rors into this distribution. For this reason, the distribution 
of observed galaxies can only give a reasonable impression 
of the shape of the LF at the bright end, which is observable 
across the entire bin. 

The observed distribution does not agree with our LFs 
at the faint end, which is not surprising, because no vol- 
ume correction has been applied; the observed distribution 
is therefore likely to be a lower bound on the true LF, be- 
cause the prevalence of faint galaxies, which are not visible 
across the entire volume of a given redshift bin, has certainly 
been underestimated. 

Nonetheless, Fig. [8] reveals a remarkably good agree- 
ment between our best-fit LFs calculated with the PML and 
the directly observed galaxy distribution, both with regard 
to the shape of the bright end and the normalization. In 
comparison, in the lowest-redshift bin, which contains the 
redshift of z — 0.1 for which the Blanton et al. LF holds, 
the observed distribution of absolute magnitudes extends to 
far brighter magnitudes than that LF would permit. 



bin 


z 




a 


M* 


lg$* 




early types 








1 


0.02 - 


0.10 


-0.639 


-21.92 


-2.502 


2 


0.10 - 


0.25 


-0.640 


-21.84t° :£ 


9 ni i +0.063 

z.on_ 063 


3 


0.25 - 


0.40 


-0.640 


-21.64±° ;fo 


9 7O9+0 063 
z ' ,u ^-0.063 


4 


0.40 - 


0.60 


-0.640 


-21.89±°;i| 


9 cqq+0.038 
i,UOJ -0.050 


5 


0.60 - 


0.80 


-0.640 


-22.0l«; 1 1 ° 


9 ^,9+0.038 
i " JU,i - 0.038 


6 


0.80 - 


1.20 


-0.640 


-21.80±8:g? 


9 7 97+0.038 
z - ' z ' -0.037 


7 


1.20 — 


1.50 


-0.639 


-21.87 


9 791 +0.051 
z -' zl -0.055 


8 


1.50 - 


2.50 


-0.639 


-21.85 


-2 694+ 071 


9 


2.50 - 


3.50 


-0.639 


-21.83 


-2.850 


10 


3.50 - 


6.00 


-0.639 


-21.80 


-2.953 




late types 








1 


0.02 - 


0.10 


-1.664 


-21.68 


-3.158 


2 


0.10 - 


0.25 


-1.637 


-20 60 +0 - 43 
zu.ou_ 50 


-2.t>ss_ 150 


3 


0.25 - 


0.40 


-1.637 


-21.62±°;p 


o n o 7 +0.100 
J ' uo ' -0.112 


1 


0.40 - 


0.60 


-1.637 


-2i.8oi° : » 


9 QQ7+0.063 
i3o '-0.075 


5 


0.60 - 


0.80 


-1.637 


-22.77t°;i^ 


n 01 9+O.O63 
J- Zlz _0 .063 


6 


0.80 - 


1.20 


-1.637 


-22.08±° f 


m 7+0.050 
°' U1 ' -0.050 


7 


1.20 — 


1.50 


-1.664 


-22.75 


Q cqq + 0.070 
o.uoo_ ogl 


8 


1.50 - 


2.50 


-1.664 


-22.97 


z. <uu_ 007 


9 


2.50 - 


3.50 


-1.664 


-22.97 


-3.252 


10 


3.50 - 


6.00 


-1.664 


-22.97 


-3.280 



Best-fit LF parameters in individual 
out in bins 2—8. Best-fit values for 
are shown for completeness, but lie 
region. 



rcdshift bins. Fits are carried 
<E>* in redshift bins 7 and 8 
beyond the reliable redshift 



The disagreement with the low-redshift LF of Blanton 
et al. is surprising. It cannot be due to a specific fault of our 
LF algorithm, because even just the absolute magnitude dis- 
tribution of spectroscopically confirmed galaxies is clearly in 
excess of the former's prediction. The fault does not appear 
to lie with our calculation of absolute magnitudes either, be- 
cause comparison to absolute magnitudes published in the 
COMBO-17 catalog shows at most small (AM « -0.06) 
systematics even as far as z ~ 1 and scatter in the recovered 
absolute magnitudes of ~ 0.22, not nearly enough to push 
the bright end forward by about a magnitude. At the ex- 
treme low-redshift end, our recovered absolute magnitudes 
are also in excellent agreement with comparison values cal- 
culated naively from the observer-frame r-band magnitudes 
simply by applying a distance modulus, indicating that our 
procedure of calculating absolute magnitudes from SED fits 
to all filter bands simultaneously does not introduce dis- 
crepancies from the conventional procedure of applying a 
fc-correction to the measured photometry in a single band. 

In this context, it i s wort h pointing out that the LFs re- 
covered by I Wolf et aT] (|2003l ) from COMBO-17 (which was 
obtained with the same telescope and instrument and over- 
laps with the E-CDFS field that we analyze here) are also 
significantly brighter than the Blanton et al. estimate and as 
bright as or even brighter than our solution. For comparison, 
we have overplotted the total LFs of Wolf et al. in Fig. [H] as 
dashed lines. In the upper left and upper right panels (cov- 
ering the redshift ranges 0.1 ^ 2 < 0.25 and 0.25 ^ z < 0.4, 
respectively, the Wolf et al. LFs shown are those at z = 0.3 
(the lowest redshift bin independently constrained by their 
work), in the bottom left panel, it is the LF at z = 0.5, 
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Figure 8. Comparison between the distribution of spectroscopi- 
cally confirmed galaxies in absolute magnitude (histogram) with 
our LF fi ts to the photometric catalog (solid line), the z=0.1 field 
LF from lBlanton et al.l J2003t) (dot ted line), and LFs in similar 
redshift bins from lWolf et al.l j2003h (dashed line). At the bright 
end, it confirms the shape of the PML solution. At the faint end, 
the observed distribution is a lower limit, because no volume cor- 
rection has been applied. 



and in the bottom right panel, z — 0.7, which coincides well 
with our choice of binning. In all bins, the Wolf et al. LF 
is substantially brighter than the estimate of Blanton et al. 
at z — 0.1, and agrees qualitatively with or is even brighter 
than our LF determination. The faint-end slopes in their 
work show considerable variance, but in the lowest-redshift 
bin, where it may be supposed to be constrained best, it 

matches ours in normalization and slope. 

The agreement between our work and that o f lWolf et al.l 

l|2003l ) with regard to the bright end at low redshift be- 
ing su bstantially brighter than estimated bv lBlanton et all 
l|2003t) . suggests that the depth of the photometry, and par- 
ticularly the surface brightness depth, may be at least partly 
responsible for th e brighter faint end recovered by us and 
IWolf et "al] l|2003l ). In particular, it is possible that the of- 
ten very extended envelopes of luminous early-type galaxies, 
which dominate the bright tip of the LF, might be underes- 
timated in the shallower SDSS data. 

Sample variance due to large scale structure must also 
be taken into account as a possible additional contributor 
to the discrepancy between our LF and that of Blanton et 
al.; Fig. [8] attests to considerable variance from one redshift 
bin to the next. We do recover a significantly fainter LF in 
the redshift range 0.25 < z < 0.4 than in the 0.1 < z < 0.25 
bin; however, this bin appears to be exceptionally faint and 
underdense. 

Variations in the high- and low-z background parame- 
ters that have so far been held constant will not shift M* to 
significantly fainter values, as we have argued in £|4.3I 



4.6 The evolution of the luminosity density 

Covariances between many of the parameters involved in 
describing the shape and normalization of the LF — such 
as between M* and S, or M* and <£>* — imply that a fairly 
large range of parameters can reproduce LFs with essentially 
similar shapes and evolutionary behaviour. For example, for 
a very steep faint end slope a, M* and $* are highly degen- 
erate, so that the same evolution with z can be described 
either as an evolution of M* or of $* . 

Specifically for the comparison to other surveys and to 
galaxy evolution models, it is thus desirable to transform 
the derived LF parameters to another quantity that does 
not suffer from these degeneracies. One such quantity is the 
luminosity density, which can be obtained by integrating the 
luminosity- weighted LF. We do so in each redshift bin and 
for early and late type objects separately. We integrate the 
LFs from Table [4] to a limit of Me, = —14. Assuming a rea- 
sonable extrapolation of the LF to fainter magnitudes, this 
limit includes almost the total luminosity density, even for 
the steep LF of late type galaxies. The results are shown in 
Fig. [9] Empty data points illustrate the luminosity density 
of late-type galaxies, and filled points of early-type galaxies. 
Error bars are calculated from the range of models within 
the la uncertainty contours in the M*-$* plane. Covari- 
ances with a are neglected, but the luminosity density is 
typically dominated by galaxies near M* , so that the actual 
contribution of the faint end to the luminosity density is 
very minor. To illustrate this, horizontal bars plotted above 
z — 0.4 show the luminosity density in each bin integrated 
only to the approximate effective absolute magnitude limit 
(determined from Fig. Q. For the early type objects, the dif- 
ference is minute up to z = 1. For late types, the difference 
is more appreciable because of the much steeper faint end; 
above z = 1, at most half of the projected luminosity density 
would be contributed by galaxies within our survey limits. 
Nonetheless, our qualitative conclusions are not affected by 
this. 

As we noted previously in t|4.4l with regard to the evo- 
lution of $* , the presence of density inhomogeneities in the 
field adds too much scatter to quantify how the luminosity 
density evolves with redshift. We will therefore address this 
question again in a subsequent analysis based on the full 
MUSYC sample, drawn from four fields in different regions 
of the sky. At z > 1.2, large jumps in the luminosity density 
of late types in redshift bins 7 and 8 reflect the problems 
that we have alluded to in H4.3I namely, that these bins fall 
outside the trusted redshift region defined in that section, 
and are only included in the optimization due to the strong 
covariances between neighbouring bins, in order not to re- 
strict the number of degrees of freedom in the search for the 
maximum likelihood solution unreasonably, while their best- 
fit values themselves are questionable. In particular, the very 
high value for the luminosity density of late-type galaxies in 
redshift bin 8 is likely to be an artifact. 

For compa rison, we have ind icated the r-band luminos- 
ity densities of I Wolf et all i|2003T ) (transformed to H = 71 
and after reclassifying their Types 1 and 2 as early types and 
Types 3 and 4 as late types) as dark (early types) and light 
(late types) shaded areas. On average, their luminosity den- 
sities are higher than ours by ~ 0.2 dex. The results agree 
very well in our redshift bin 5, which comprises one of the 
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Figure 9. Luminosity density as a function of rcdshift for early 
types (filled diamonds) a nd late type s (emp ty diamonds). Shaded 
regions show values from I Wolf et al.l 1120031 ) for comparison (light 
for late, dark for early types). Data points above z=1.2 are in- 
cluded in the optimization, but lie beyond our pre-defined rcdshift 
range of confidence. Horizontal bars show the luminosity density 
integrated approximately to the effective magnitude limit. 



known large scale structure features, but their data do not 
show the strong underdensity in the region 0.25 ^ z < 0.4 
as pronouncedly as ours do. 



4.7 The detailed shape of the LF 

A striking feature of our LF solution is the extremely steep 
faint end of the LF of late- type galaxies of a « —1.62. 
This steep slope dominates the total LF at the faint end, 
leading to an upturn at intermediate magnitudes. How- 
ever, even this double-Schecbter shape may not provide a 
comprehensive description of the faint end of the LF, be- 
cause even within the subclasses of early and late types 
that we analyze, features that are not consistent with the 
Schechter functions, such as up- or d ownturns, may exist . 
In particular, galaxy evolution mode ls (|Springel et alJfeoOlt 
|Pe Lucia. KaufFmann fc Whitell2004h . which propose a host 
of mechanisms for suppressing star formation in low mass 
halos, and the relatively low number counts of ultra-faint 
dwarf galaxies in the Local G roup (SPSS), which show good 
agreement with such models dMaccio. Kang fc Mo"orell2009l ; 
IKoposov et~aI1l200g| ). may require a break in the late-type 
LF if they are to be brought into consistency with our claim 
of a very steep faint end slope at intermediate magnitudes. 
This raises two questions: Down to what absolute magni- 
tude can our analysis constrain the LF, and is our best- 
fit Schechter function really an adequate fit over this entire 
range, or are there indications of a departure from this shape 
at faint luminosities? 

This question is answered best by representing our LF 
by a binned distribution, rather than by a low-parameter an- 



alytical ansatz. As in §4T3] we solve for the LF over the range 
0.1 ^ z < 1.2, but instead of modeling the SED-specific LF 
with evolving Schechter functions, we describe it with se- 
ries of bins in absolute magnitude of 1 mag width each. The 
bin boundaries are taken relative to M* in a given redshift 
bin, and the normalization of the LF value <& in each mag- 
nitude bin is also relative to the normalization parameter 
$* in a redshift bin, as inferred from our solution in Ta- 
ble [3] Therefore, the binned solution represents the shape of 
the LF after removing the evolution in M* and <&*. All LF 
parameters not explicitly represented by this binned ansatz 
are held constant at their best-fit values from Table 

For our first analysis, we scan the likelihood distribu- 
tion for one magnitude bin at a time as a function of $ 
(the value of the LF across this bin) , while keeping all other 
bins at the values implied by the Schechter function. The 
resulting uncertainty range is a measure of how much an 
individual magnitude bin may deviate from the Schechter 
fit, and will thus provide very conservative upper and lower 
limits on the shape of the LF at the faint end. If these limits 
are sufficiently tight, they will allow us to rule out major- 
departures from the Schechter function. 

Fig. [10] shows the results of this calculation. The LF 
has been adjusted to the M* and $* values in redshift bin 
2 from Table Thin and thick solid lines show the best- 
fit Schechter functions, as calculated previously, for late and 
early types, respectively, and the dotted line illustrates their 
sum, which is dominated by early times at the bright end 
and by late types at the faint end. Data points (filled for 
early types, empty for late types) show the maximum likeli- 
hood value for each bin, and error bars the acceptable 68% 
confidence range for one degree of freedom. 

The superposition of the shallow early and steep late 
type LFs produces a distinct upturn at intermediate magni- 
tudes of M r sa —19. Such an upturn has been suggested by 
iBlanton etafl (|2005l ") on the basis of SDSS data. In the latter 
survey, the fairly shallow surface brightness limit prevents 
a direct measurement of the full faint end of the LF, but 
extrapolating for galaxies beyond this limit, their steepest 
estimate for the faint-end component is a = —1.52, slightly 
shallower than our result. This estimate is shown in Fig. 1101 
by a dashed line. Around the onset of the faint-end compo- 
nent, at M r « —19, our LF agrees very well in normaliza- 
tion and shape with the LF of Blanton et al. The Blanton 
LF does deviate significantly from ours at the bright end, a 
phenomenon already discussed in £|4.5I 

The principal conclusion from this figure is that we can 
constrain the late- type LF to M r ~ —14, and the early- type 
LF to M r ~ —15. A downturn in either LF above these 
magnitude limits can be ruled out. If mechanisms suppress- 
ing star formation do act to lower the abundance of dwarf 
galaxies from the extrapolated a » —1.62 LF, this must 
happen at even fainter absolute magnitudes. There is some 
room for an excess of faint galaxies over the Schechter fits. 

To gain a more accurate picture of the shape of the LF, 
we now carry out a simultaneous optimization of all magni- 
tude bins of both SEDs. After finding the best-fit solution, 
we obtain uncertainties with the previously-described strat- 
egy (for each free parameter, we scan the range of values 
that yield an acceptable offset in \ 2 from the best-fit solu- 
tion, allowing all other magnitude bins to be optimized in 
the background). The resulting LF is shown in Fig. 1111 Er- 
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Figure 10. Binned representation of the shape of the late- and 
early-type LFs at 0.1 $C z < 1.2, shifted to M* and <j>* in red- 
shift bin 2. Magnitude bins have been optimized one at a time. 
The thick solid curve represents the best-fit Schechter function of 
early types, the thin curve of late types. Data points with error 
bars show the uncertainties on $ in any one individual magnitude 
bin under the (idealized) assumption that all bins are indepen- 
dent. They serve as an upper limit on any deviation of the LF 
shape from the Schechter function, and show that the late-type 
LF does not exhibit a break down to at least M r fs —14, but 
instead continues to rise steeply with a —1.65. The dotted line 
represents the sum of both LFs, and th e dashed line the steep- 
est LF suggested bv lBlanton et al I d2005h . which agrees well with 
ours around the onset of the faint end upturn, although it fails to 
match our LF at the bright end. 

ror bars are again calculated for one degree of freedom and 
a confidence level of 68%, as appropriate for an error calcu- 
lation in a single parameter. The uncertainties are generally 
smaller than in the previous analysis, which at first appears 
counterintuitive (allowing the "uninteresting" parameters to 
find a conditional maximum likelihood generally enlarges the 
uncertainties on any one "interesting" parameter). However, 
this can be explained by the fact that we are now evaluating 
the uncertainties relative to the true maximum likelihood fit 
in the parameter space consisting of the individual magni- 
tude bins, rather than the best Schechter fit. 

This binned LF is sensitive to our assumptions for the 
values of the constant "background" LF parameters at very 
low and very high redshift. To illustrate this sensitivity, we 
have taken the output from the Monte Carlo trials discussed 
in 14.31 and calculated the binned LF for those trials not re- 
jected due to their low likelihood. The shaded areas in Fig. 
[TT]show the full range covered by these binned LFs. 

The figure confirms the basic picture from Fig. 1101 The 
late- type LF continues to rise steeply to M r w — 14 with- 
out any indication of a break. Of all our Monte Carlo trial 
runs with acceptable likelihoods, not one exhibits a strong 
downturn in the faint end slope at M r < —14. Although 
the last two significant bins in the late-type LF do exhibit 




Figure 11. Binned representation of the shape of the late- and 
early-type LFs at 0.1 $J z < 1.2, shifted to M* and <!>* in redshift 
bin 2. Here, all magnitude bins have been optimized simultane- 
ously, representing our best estimate for the shape of the LF with 
associated uncertainties around the best fit. Shaded areas show 
how variations in the high- and low-redshift background LF pa- 
rameters propagate into the best-fit solution. 

small, but significant deviations from the Schechter fit, the 
Monte Carlo results show that variations in the high- and 
low-redshift background can account for their magnitude. 

The early-type LF is constrained well down to M r = 
— 16, and upper limits can be applied down to M r = —14. 
These show that the LF follows our Schechter fit to at least 
M r = —16. The uncertainties leave room for deviations from 
the power-law faint end slope in either sense at fainter mag- 
nitudes, but nonetheless the late-type LF clearly dominates 
the total LF at least to M r < -14. 

We have repeated this analysis using the more strin- 
gent star/galaxy classifier limit of 0.8, instead of our default 
0.97. This incurs the cost of possibly excluding very com- 
pact dwarf galaxies, but is closer to typical values adopted 
in the literature. As we have seen in 34.31 the best-fit early- 
type LF is characterized by a slightly shallower faint-end 
slope and fainter M* , while the late-type LF is almost 
unaffected. However, calculation of the binned LF in the 
manner described above shows no qualitative difference in 
the shape of the LF, compared to our default calculation, 
which indicates that the recovered LF shape does not de- 
pend on the precise star/galaxy threshold (at least in the 
range 0.8 ^ S/G ^ 0.97, which is spanned by these two 
samples) . 



5 CONCLUSIONS 

We have proposed a new Bayesian method for the calcu- 
lation of galaxy luminosity functions from multi-band pho- 
tometric samples. The method, the Photometric Maximum 
Likelihood Method, deconvolves the observed distribution of 
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galaxies as a function of fluxes in multiple filter bands into 
the luminosity functions of all galaxy populations (identi- 
fied by redshift and SED type) that may contribute to this 
observed distribution. It thus eliminates the requirement to 
assign photometric redshifts to individual galaxies. This ap- 
proach offers multiple advantages: 

• The PML achieves a proper deconvolution of the ob- 
served data into the constituent galaxy populations. 

• The entire information contained in the multi-band 
photometric data is retained for the calculation of the LF, 
instead of incurring the information loss that is associated 
with reducing these photometric data to a single best-fit 
photometric redshift. 

• There is no requirement to model the often very com- 
plicated error functions of photometric redshifts; instead, 
only an understanding of the photometric error function is 
required. 

• Conventional photo-z determinations and LF calcula- 
tions based on them are dependent on many (implicit or 
explicit) priors about the galaxy distribution as a function 
of redshift, which often go unacknowledged. By solving for 
the parameters of multiple galaxy populations over a large 
range of redshifts simultaneously, the PML significantly re- 
duces the number of priors and allows for a self-consistent 
solution. 

We have applied this algorithm to the E-CDFS field 
of MUSYC (the Multi- Wavelength Survey Yale-Chile) with 
the aim of recovering the faint-end shape of the LF at low 
redshift. For this purpose, we have represented the galaxy 
population up to z = 1.2 by two SED-specific Schechter 
functions (one for early-, one for late-type galaxies). Princi- 
pal results are: 

• In accordance with the current consensus in the litera- 
ture, the LF of early-type galaxies is characterized by a fairly 
bright characteristic absolute magnitude M* and a shallow 
faint end slope, while the LF of late- type galaxies is much 
steeper. 

• The faint-end slope of the late-type galaxy population 
may be as steep as a ~ —1.62. The superposition of these 
two SED-specific LFs produces an overall LF that cannot 
be described by a single Schechter function, but instead ex- 
hibits an upturn at intermediate magnitudes, where the late- 
type LF begins to dominate the overall LF. This upturn is 
noticeable in comparison to a standard a — —1.2 LF at 
magnitudes fainter than M « M* + 3 (M w —19 in our 
sample) . 

• The bright-end exponentia l cutoff is significantly 
brighter than found in the SDSS (iBlanton et al.ll2003h. bu t 
slightly fainter than found in COMBO-17 (|Wolf et aljl2003h 
from photometric redshifts. 

• We find fluctuations in the normalization of the LF 
and in the luminosity density with redshift in the E-CDFS 
that are broadly consistent with large-scale density inhomo- 
geneities known from spectroscopic data but prohibit firm 
conclusions regarding any systematic trend of these quanti- 
ties with redshift. 

• Calculation of the shape of the LF using a binned 
ansatz, rather than the analytical Schechter form, shows that 
the steep faint end extends at least as faint as M r ~ —14, 
which confirms the presence of a faint end upturn. Therefore, 



we postulate a significantly larger number of dwarf galaxies 
in the universe than would be predicted by simply extrapo- 
lating the typical faint end slope of a ps —1.2 recovered in 
shallower field galaxy surveys to fainter magnitudes. 

• Our results show that, although our conclusions are 
qualitatively robust, statistical results still depend on pri- 
ors, such as the extreme high-redshift LF, that are difficult 
to constrain in a self-consistent way, even with a high-quality 
data set such as MUSYC. Any work attempting to constrain 
these quantities from purely photometric data must give a 
careful account of the effect of such priors on its conclusions. 

The PML formalism offers many opportunities for re- 
finements. It can be easily generalized to calculate distri- 
bution functions over colors, stellar masses, or star forma- 
tion rates, as long as these quantities can be unambiguously 
assigned to the SED templates used in the analysis. To im- 
prove the calculation of LFs at high redshift, it is possible to 
implement passive redshift evolution in the SED templates, 
as opposed to the redshift-zero empirical templates which 
we have used in our analysis. We are planning a subsequent 
study of the full MUSYC catalog which will explore some of 
these avenues. 
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APPENDIX A: MODIFICATIONS TO THE 
PUBLISHED MUSYC PHOTOMETRIC 
CATALOG 

We have made two modifications to the published catalogs 
of the MUSYC-ECDFS in order to optimize the accuracy 
of the photometric redshifts, and, by extension, of our LF 
results. The first of these is a zero-point correction. Zero- 
point offsets manifest themselves in photometric fluxes in 
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one given filter band generally being determined higher or 
lower than the true fluxes by a constant factor, and can be 
due to an inaccuracy in the photometric calibration or the 
transmission curves. In order to detect such offsets, we have 
carried out SED fits to all galaxies with known spectroscopic 
redshifts. The residuals of the observed fluxes relative to the 
fluxes predicted from the best-fit SED are an approximation 
of the zero-point offset of the fluxes recorded through this 
filter. By visual inspection of the distribution of photometric 
residuals, we have determined the required corrections to 
the published fluxes as a function of the filter band and list 
them in Table [ATI The fluxes are corrected according to the 
formula 

fcorr = /published * 10 • (Al) 

This procedure is also sensitive to template insuffi- 
ciency, i.e., an inability of the template set to reproduce 
the entire range of SEDs displayed by our sample objects. 
Ascribing systematic photometric residuals that are really 
caused by template insufficiencies to photometric zero point 
shifts instead may optimize the quality of photo-z fits to the 
spectroscopic subsample, but possibly introduce systematic 
errors for other galaxy populations and other redshifts. In or- 
der to provide an independent verification of our zero-point 
shifts, we have applied a second test that is less sensitive 
to template insufficiencies: We have compared photometric 
residuals between two different filters, but for the same SED 
and rest-frame wavelength, i.e., we have examined the dis- 
tribution of the quantity 

SA = Ig |i -Ig f^- | 2 = lg A(filten) - Ig A{filter 2 ) 

Jpred Jpred 

+lgc 2 (SED,\r est ) [i -lgc 2 {SED,\r es t) | 2 

for all possible pairings of galaxies fitted with the same 
SED template and as a function of the difference in rest- 
frame wavelength. The term C2 is a hypothetical correction 
to the SED template, dependent on the template itself and 
the rest-frame wavelength. If the difference in A res t between 
the two observations is small and both are drawn from the 
same SED, the third and fourth logarithms cancel, and this 
becomes a measure of the difference in photometric resid- 
uals, evaluated for the same SED type at the same A rest , 
but evaluated through different filters at different redshifts. 
In practice, the corrections C2 will depend on the true SED 
of the observed objects, which may deviate from the best- 
fit SED template, so that the logarithms that depend on 
C2 do not cancel perfectly, but the size of such deviations 
is limited by the fact that both objects are fitted by the 
same SED template, so that this method allows for at least 
a qualitative verification of the previously found zero-point 
shifts. 

The resulting relative zero-point shifts between any two 
optical filter bands are consistent with the numbers pre- 
sented in Table |A"T1 bolstering our confidence that the pho- 
tometric residuals can indeed be predominantly ascribed to 
zero-point adjustments, rather than template insufficiencies. 
We do point out that we have not been able to use this 
method to independently verify our zero-point shifts for the 
J- and A'-bands, because the number of objects permitting 
us to observe the same Arest through both filters at different 
redshifts is too small. However, the actual impact of these 



Table Al. Adjustments to Photometric Catalog 



Filter 


A [mag] 


ei 


e2 


U 


-0.075 


0.1 


0. 


B 


+0.00 


0.05 


0. 


V 


+0.0125 


0.04 





R 


+0.025 


0.03 


0. 


i 


-0.1 


0.08 


0.1 


z 


+0.075 


0.09 


-0.4 


.] 


-0.15 


0.15 





K 


-0.075 


0.2 


0. 



Parameters for corrections to the published photometric fluxes 
and uncertainties, according to Eg. lAll and Eg. IA3I 

zero-point corrections on our LF results is minor and does 
not alter our qualitative conclusions. 

The second adjustment to the photometric catalogue 
concerns the catalog uncertainties. We have tested the cat- 
alog MUSYC photometric uncertainties by comparing the 
scatter of the residual photometric fluxes from SED fitting 
to objects with known redshifts. While, in the ideal case, 
the residuals should show a Gaussian distribution (i.e., 68% 
should fall within the published standard deviation, 95% 
should fall within two standard deviations, etc.), the residu- 
als in the MUSYC catalog typically exhibit a slightly larger 
scatter, indicating that the catalog uncertainties are an un- 
derestimate, and the scatter increases with increasing flux. 
The latter behaviour is easily understood, because the pub- 
lished uncertainties only reflect internal photometric consis- 
tency, but do not include photometric calibration uncertain- 
ties or local background errors. 

In order to achieve a more realistic, Gaussian distribu- 
tion of the uncertainties, we have introduced two additional 
parameters to transform the published errors: 

cr„ew = V (cFpubiished x (1 + e 2 )) 2 + (ei x flux) 2 (A3) 

The effect of these two parameters ei and e 2 is to in- 
crease the uncertainties if they have been underestimated 
and to add a component to the uncertainty that is propor- 
tional to the flux (i.e., a zero-point calibration error). We 
have adjusted these two parameters manually for each filter 
band in order to achieve a more Gaussian-like distribution 
of the photometric residuals. The values are listed in Table 
[AH 



